NASA TECHNICAL nasa tm x- 62,309 

MEMORANDUM 


8 

CO 

s' 


< 

CO 

< 


(NASA-TH-I-62309) GEOHETBI PROGRAM FOE N74-11810 

AERODYNAMIC LIFTING SURFACE THEORY (NASA) 

108 p HC $7,50 CSCL 01A 

Unclas 

G3/01 23014 


GEOMETRY PROGRAM FOR AERODYNAMIC 
LIFTING SURFACE THEORY 

Richard T. Medan 


Ames Research Center 
Moffett Field, Calif. 94035 



September 1973 



TAELE OF CONTENTS 

1. INTRODUCT7G N 

2. NOTATION 

2.1 ARRAY NOTATION 

3. INPUT AND OUTPUT DESCRIPTION 

3.1 INPUT DATA 

3.2 GEC MET BY PILE 

3.3 BOUNDARY CONDITION FILE 

4. METHOD OF SOLUTION 

5. USE OF TRANSFORMATION CAP ABILITY 

6. SAMPLE CASES 

7. CONTROL CARDS 

7.1 AMIS' IBM T3S SYSTEM 

8. REFERENCES 
APPENDIX A 
FIGURES 

OUTPUT FOR TEST CASES 
PROGRAM LISTING 


i 




GEOMETRY PROGRAM FOR AERODYNAMIC 


LIFTING SURFACE THEORY 
Richard T, Medan 
Ames Research Center 


ABSTRACT 


This document is a description of and user's manual for a USA FORTRAN 
IV computer program that provides the geometry and boundary conditions 
appropriate for an analysis of a lifting, thin wing with control surfaces in 
linearized, subsonic, steady flow using kernel function method lifting surface 
theories. The data which is generated by the program is stored on disk files 
or tapes for later use by programs which calculate an influence matrix, plot 
the wing planform, and evaluate the loads on the wing. In addition to pro- 
cessing data for subsequent use in a lifting surface analysis, the program is 
useful for computing area and mean geometric chords of the wing and control 
surfaces. Any planform shapes, including asymmetrical ones, ones with 
mixed straight and curved edges, and ones with control surfaces may be 
handled. The program is able to compute a control surface downwash mode 
and residual mode and a control surface lifting pressure mode. This in turn 
allows conventional lifting surface methods to effectively handle partial span, 
trailing edge control surfaces. 
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GEOMETRY PROGRAM FOR AERODYNAMIC 
LIFTING SURFACE THEORY 


Richard T . Medan 
Ames Research Center 
September 1973 


1 I NTR OD UCTI ON 


This document is a description of and user's manual for 
a USA FORTRAN IV computer program that provides the 
geometry and boundary conditions appropriate for an 
analysis of a lifting, thin wing with control surfaces 
in linearized, subsonic, steady flow using kernel 
function method lifting surface theories. The data 
which is generated by the program is stored on disk 
files or tapes for later use by programs which 
calculate an influence matrix, plot the wing planform, 
and evaluate the loads on the wing. In addition to 
processing data for subsequent use in a lifting surface 
analysis, the program is useful for computing area and 
mean geometric chords of the wing and control surfaces. 
Any planform shapes, including asymmetrical ones, ones 
with mixed straight and curved edges, and ones with 
control surfaces may be handled. The program is able 
to compute a control surface downvash mode and residual 
mode and a control surface lifting pressure mode. This 
in turn allows conventional lifting surface methods to 
effectively handle partial span, trailing edge control 
surfaces (ref. 1) . 

Questions concerning either this document or the 
computer program or the associated computer programs 
should be directed to 

R. T. Medan 
Mail Stop 221-2 
Ames Research Center 
Moffett Field, 

Ca. 94035 
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•GE CHET BY PROGRAM 


2 NOT ATTCN 


2.1 ARRAY NOTATION 

This notation will be explained by example. 

(ETA) ®eans ETA is an array 

(ETA(JJMAX)) refers to the elements of (ETA) 

from 1 through JJHAX 

(ALFAB (F P,NF) ) refers to the first P? rows and 

NF columns of the two 
dimensional array (ALFAB) 

ALFAB(3,NF) refers to the element of (ALFAB) 

in row 3 and column NF 
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GEOMETRY PROGRAM 


3 IN PDT AND O UTP UT DES CRIPTION 


The function of the program is to determine and store 
on data files certain geometry information and boundary 
conditions for a wing planform given a very basic 
description of the planform. The specific information 
to be determined and stored on the data files, which 
are used by other computer programs, is discussed in 
sections 3.2 and 3.3 fallowing a description of the 
input in section 3.1. The coordinate system is defined 
in fig. 1. 


3. 1 . INPUT DATA 


The input data is described and defined in the 
following paragraphs in the order in which it is 
read by the program. The input FOBMAT is 8F10.O for 
floating point numbers, 1615 for integer numbers, 
and 10L1 for logical variables unless otherwise 
indicated. Each item in the following list 
represents one or more cards and, depending on the 
data of previous items, may or may not be read by 
the program. All coordinates and areas can be given 
in any consistent units. 


3.1.1 ITEM 1 
JTITIEl 

This array is used for titling information. The 
information may consist of up to 80 characters 
(somewhat less if the integer word length of the 
computer does not divide evenly into 80) . if the 
planform plotting program is to be used, the 
title should be centered in the first 42 columns. 
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GEOMETRY PROGRAM 


3. 1.2 ITEM 2 

ID 

This is an Identification number which can either 
be given by the user or assigned by the program. 
If the input value is zero, the program will read 
an identification number file, increment the 
first value read by one, rewrite the 
identification number file using the incremented 
value, and use the incremented value as the 
identification number. In order to do this a 
file must already exist with four integer 
identification numbers on it in unformatted 
(binary) form. The unit numher of this file is 
determined by the value of the inteqer »9 in the 
program (see the DATA statements) . 

C WTYP E (cols. 6-10) 

This is an integer denoting the type of chordwise 
control point distribution for applying the 
boundary conditions. 

<0 indicates that the program should accept a 
control point distribution given by the 
user. See section 3.1.21. 

=0 indicates that the program should put 
ccntrol points first on the trailing edge, 
then on the leading edge, then 
equally— distributed by the cosine rule on 
the interior (fig. 2). 

>0 indicates that the program should compute 
and store ccntrol points using a modified 
cosine rule (fig. 3). 

SWTYPE (cols. 11-15) 

This is an integer denoting the type of spanwise 
ccntrol point distribution. 

=0 indicates that the program should compute 
and store a uniform cosine distribution 
(fig. 4) • 

indicates that the program should accept and 
store a user-defined distribution * (see 
3.1.22). 
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GECMETBY PRCGBAB 


IlllfC (cols. 16-20) 

A non— zero value causes the integration stations 
and information at these stations to be printed 
(see test cases) . 

C PRES S (cols. 21-25) 

A non— zero value causes the flap lifting pressure 
mode (provided other input indicates that a flap 
exists (see 3.1.3) ) to bo computed at NC (see 

3.1.23) chordwise stations for each spanwise 
ccntrol station. 

DNSYM (cels. 26-30) 

This is an integer. It should be zero if and 

only if the leading and trailing edges are 
symmetric with respect to the center line. 

TRNFRM (cols. 31-35) 

This is an integer indicating the type of 

transformation to be performed upon the wing edge 
definitions input by the user. The 

transformation will only be dene if TR8ANG (see 
3.1.6) is not zero. See Chapter 5 for 

considerations relevant to the use of the 

transformation capability. 

<0 indicates a cranking transformation (fig. 5) 

=0 indicates a shearing transformation (fig. 5) 

>0 indicates a rotating transformation (fig. 5) 

CUT (cols. 36-40) 

This is an integer indicating the output, level 
for SUBROUTINE INTGRT, the program which computes 
wing and flap areas and mean geometric chords. 

-1 gives nothing at all 

0 gives warning messages 

1 gives the above + answers 

2 gives the above + the given and effective 
locations of planform kinks 

3 gives the above + input function values and 
integration stations 

N07 (cols. 41-45) 

A non- zero value suppresses writing the geometry 
(section 3.2). A zero value causes it to be 
wr iften on unit W7 (see the DATA statements). 
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GEOMETRY PROGRAM 


N08 {cols. 46-50) 

A non-zero value suppresses writing the boundary 
condition file (section 3-2) • A zero value causes 
it to be written on unit S8 (see the DATA 
statements) . 


3.1.3 ITEM 3 

JBCSJJOU 

This is a logical array specifying which 
symmetric boundary conditions will be computed 
and/or read and stored on unit H8 (see 3.3). A 
full explanation of this item and also of the 
following item is not included because the 

boundary condition program which was developed 
after this one performs this function and is more 
convenient to use than the geometry program. See 
the boundary condition program listing for an 
explanation of the significance of each item of 
(BCS) and {8CAS) . BCS { 1) = . TRUE. (i.e. a T in 
column 1) refers to the uniform dcnvwash case, 
causing this case to be computed and stored 
(NC8=0) . BCS (2) =- TRUE. refers to the pitching 
mode from which guasi-steady pitching derivatives 
may be obtained. BCS (5) =. TRUE. refers to the 
flap downwash mode and tells the program to 
expect flap data. 


3.1.4 ITEM 4 

(BCASliOll 

This is a logical array specifying which 
anti-symmetric boundary conditions will be 
computed and/or read and stored on unit H8 (see 
3.3). BCAS (2) =. TRUE. refers to the rolling mode 
from which quasi-steady rolling derivatives may 
be obtained. 
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GECBETBY PROGRAM 


3.1.5 ITEM 5 
BAC H 

This is the Mach number tor the flap pressure 
mode calculation. It is also the default Mach 
number for the influence matrix program and will 
be used if no other Mach number is specified at 
the time of running the influence matrix program. 

3.1.6 ITEM 6 (FORMAT = 2I5,3F10.0) 

KL 

This is the number of pairs of points on the 
entire wing for definition of the leading edge. 
The maximum value allowed is 100. 

IETYPE {cols. 6-10} 

This variable controls the type of interpolation 
to be used hetween the data points for the 
leading edge. 

=0 Linear interpolation between points. 

TO A controlled deviation interpolation will be 
used between points (ref. 2). This method 
can handle edges with mixed straight and 
curved sections. To input a section which 
is a straight line, give 3 points lying on 
the line. 

LESCAL (cols. 11-20) 

This is a scale factor (floating point) for the 
XL’s. The default value is 1.0. 

X£M (cols. 21-30) 

This is an interpolation control constant used 
when LETYPRT0. A value of 0. will give linear 
interpolation in the end intervals. A value of 
1.0 will give parabolic interpolation in the end 
intervals. A value in between will give a curve 
in between. The end intervals are the spaces 
closest to the wing tips. 

TRNANG (cols. 31-40) 

This is the cranking, shearing, or rotation angle 
to be imposed on the planform data. The angle 
should be in degrees. Input 2 ero if nc 
transformation is desired. 
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GEOMETRY PROGBAH 


3.1. 7 ITEM 7 

(YL ( NMA X) l 

This”array is to contain the spanwise coordinates 
of the leading edge definition points. Only the 
points for Y>0 need be given when the leading 
edge is symmetric with respect to Y. The values 
must be in descending crder and Yl(1) must be the 
spanwise coordinate of the right tip. If 
UNS YM-0 , then N M A X = (NL + 1) / 2 . If USSYMtO, then 
NMAX=NI. 


3.1.6 ITEM 8 

(XI (FMA X) ) 

This array is to contain the streamwise 
coordinates of the leading edge definition 
points. These values will be multiplied by the 
scale factor LESCAL immediat ely after being read 
in. For the definition of Nil AX see ITEM 7. 


3.1.9 ITEM 9 — read only if LETYPEf 0 
NDL 

The number of Y stations of ITEM 8 cn the entire 
wing where the leading edge is kinked. The 
program needs to know where the kinks are in 
order to integrate for the wing area and mean 
geometric chord mere accurately. 


3.1.10 ITEM 10— read only if LETYPEfO and NDLfO, 

(INDDL ( RMAX ) ) 

This is to be an integer array giving the Y 
stations where there is a kink. A value of 2 
would, for example, mean that the second number 
of ITEM 7 would be the location of a kink. If 
the wing has a symmetric effective leading edge, 
cnly the indices corresponding to Y>0 are 
required. If UNSYHfO, then NMAX=NDL. if 
UNSYM=0, then NMAX = (NDL+1)/2. 
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GECMBTBY PRCGFAM 


3.1.11 IT5PI 11 (FORMAT- 215, 2F10.0) 

NT 

This is the number of pairs of points on the 
entire wing for definition of the trailing edge. 
The maximum value allowed is 100. 

TETYPE (cols. 6-10) 

This integer variable controls the type of 
interpolation to be used between the data points 
for the trailing edge. It corresponds to LETYPE 
cf ITEM 6. 

1ESCAL (cols. 1.1-20) 

This is a scale factor for the XT's. The default 
value is IESCAL, not 1.0. 

XCTE (cols. 21-30) 

This is the interpolation control constant for 
the trailing edge when TETYPEjO. It corresponds 
to XCLE of ITEM 6. 


3.1.12 ITEM 12 

This is the trailinq edge counterpart of ITEM 7. 
(YT) must be in descending order. YL{1) defines 
the wing tip, so YT(1) could be greater or less 
than YL ( 1 ) • N«AX=NT or NM AX- (NT + 1) /2 according 
to the value of ONSYM. 


3.1.13 ITEM 13 

JIIJ182IU 

This is the trailing edge counterpart of ITEM 8. 
The values are scaled by TESCAI just after being 
read in. 
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3. 1.14 


3.1.15 


3.1.16 


GECHJ5TRY PROGRAM 


ITEM 14 — read only if TETYPE^O 
NDT 

This is the trailing edge counterpart of ITEM 9. 


ITEM 15 
(INDDT (N M A X) ) 

This”is the trailing edge counterpart of ITEM 10. 
NM&X = NDT or NMAX= (NDT + 1) /2 according to the value 
Of OKSYR. 


ITEM 16 
AWIjJG 

This is the wing reference area and is optional. 
If not given, the program will compute the actual 
area. 

CBASW (cols. 11-20) 

This is the wing reference chord and is optional. 
If not given, the program will compute a 
reference chord. The value computed is the 
integral of c**2*dY over the span divided by the 
reference area. 

ft FI A E (cols. 21-30) 

This is the flap reference area. It also will be 
computed if not qiven and if there is a flap (see 

3.1.3) . 


CEASF (eels. 31-40) 

This~is the flap reference chord. If not given, 
and if there is a flap, then it will be computed 
as the integral cf cf**2*dY over the flap span 
divided by the flap reference area. 
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GECHETBT PH CGB A E 


3.1.17 ITEM 17 — read only if there is a flap. 

YFJ 

This is the spanwise position of the left flap 
edge. 

YF2 (cols. 11-20) 

The spanwise position cf the right flap edge. 
IF2>YF 1 . 


3*1.18 ITEM 18 — read only if there is a flap. 

CHI J 

The chord wise position of either the left flap 
corner or the hingeline at Y=0 relative to the 
local chord and nondimensionalized by the local 
chord. 0<CHI1<1« If the wing is uusynmetric 
after being transformed, then Cfill applies at. the 
left flap corner. Also, if YF1*-YF2 (indicating 
that there are two flaps) , then CHIl applies at 
the left flap corner. Ctberwise there is only a 
single flap and CflII applies at the centerline, 
while CHI2 (below) applies at both corners. In 
the latter instance, the hingeline may be kinked 
at the centerline. 

CBI2 (cols. 11-20) 

The chordwise position of tho right flap corner. 


3.1.19 ITEM 19 
JJH AX 

This is the maximum number cf spanwise 
integration stations. This number will be the 
upper limit for nest of the integrations involved 
with the current wing. It will be used to 
determine the wing and flap areas and chords. It 
will be used for the maximum number of spanwise 
integration stations in the influence matrix 
program. It will be used to integrate for the 
pitching moment in the forces program. And it 
will be used elsewhere. The value should be 
chosen as large as possible. If it is input as 
2Gro » then the program will determine the largest 
possible value net exceeding JJHXMX (see TATA 
statements) such that the following relation will 
be satisfied: JJMAX+ 1= (NN + 1) *2**K, where K is 
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3.1.20 


GEOPlBTEY PFOGEAM 


same integer. NN 
relationship between 
the control points 
integraticn points. 


ITEM 20 
FP 

This is the desired number of chordwise control 
points. The chordwise control points determined 
by the geometry program serve only as a default 
set. They may be changed at will in the 
influence matrix program. The maximum value for 
PP is 15. 

NF (cols. 6-10} 

This is the desired number of spanwise control 
points. The spanwise control points determined 
by the geometry program serve only as a default 
set. They may be changed in the influence matrix 
program, but net completely arbitrarily as in the 
case of the chordwise control points. In the 
influence matrix program (and in this prooram) 
the spanwise control points must be a subset of 
the spanwise integration points. The maximum 
value for NF is 47. 

NN (cols. 11-15} 

This is a reference number for computing the 
spanwise control points. The default value is 
NF. See 3.1.2 for further explanation. 

XL (cols. 16-20) This is the number of times 
that the spanwise integrations will be done for 
the flap pressure mode. The default value is 4. 
♦♦NOTE: currently the flap pressure mode 
calculation program is under revision and this 
capability is not available. Also it is more 
convenient to use the boundary condition program 
for this purpose***^ 


is given below. This 
NN and JJMAX occurs because 
must be a subset of the 
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GECHETBY PBOGBAM 


3.1.21 ITEM 21— read only if CWTYPECO. 

JCmCPiPPil 

Values of the chcrdwise control point, location 
referenced to the local leading edge, an 
ncndimensicnalized by the local chord. 
(KCHICP |P) <1 . Values must be in ascending order. 


3.1.22 ITEM 22 — read only if SWTYPE=fO. 

(HI ND EX (N M AX) ) 

This array is to contain integers for computing 
the spanwise control point positions by 
ETA (N) =COS (NINDEX (N) *PI/ (NN+1) ) where *1=3.14... 
If SWT YPE=0 , then the program will compute 
(NINDEX) as 1,2, 3,4,.... N M A X = N F for an 
unsymmetric wing and NM AX= (NF + 1 ) /2 for a 
symmetric wing 


3.1.23 ITEM 23 — read only if there is a flap. 

SC The number cf chordwise points at which the 
flap net pressures will he computed. The 
spanwise stations that the program will use are 
the spanwise control stations. The default value 
for NC is 49. The maximum allowable value is SO. 


3.1.24 ITEM 24 — read only if there are some downwash 

distributions to be read from cards. 
This item is read for N = 1, NMAX. If 
PP<8, then NMAX cards are read. If 
9<PP<16, then 2*NMAX cards are read. 
UMAX = (EF+1)/2 for symmetric cases on a 
symmetric wing. N M AX - NF/2 for 
antisymmetric cases cn a symmetric wing. 
UMAX = NF for an unsymmetric wing. 
Input from right tip towards left tip. 

(A1FA B (P.N1 . P= 1 . FF ) 

Downwash values~to be stored on the boundary 
condition file. Input from leading edge towards 
trialing edge. The user will have tc examine the 
program in order to use this capability of 
supplying downwash modes. The boundary condition 
program performs this function much better 
anyway. 
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GEOMETRY PROGRAM 


3.1.25 ITFM 25 
CEL T AO 

The chordwise integration accuracy parameter. 
This card is read by SUBROUTINE FLPPSiN and may be 
repeated in order to assess the convergence. Use 
the smallest value last. A value of zero causes 
a return to the main program. DELTAO=the square 
of the aspect ratio is the default value. 
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3-2 GEOHETBY FILE 


The following data is written on the file referenced 
by the integer W7 in the program. The data is 
written in binary form and comprises 2 records if 
there is no flap and 3 records if there is a flap. 
The writing of the data on the file can be 
suppressed by inputting NC7 as non-zero (see 3.1.2). 
All lengths written on this file have been 
ncndimensicnalized (usually by the effective 
semispan, 52) . 


3.2.1 FIRST LOGICAL BECCBD 

ID 

Identification number. See 3.1.2. 

FP 

The number of chordwise contrcl points. 

JNF 

The number of spanwise control points. In 
subsequent programs this number is given the name 
EH. 


CBTY PE 

Integer denoting the type of chordwise control 
point distribution. See 3.1.2. 


SBTYFE 

Integer denoting the type cf spanwise control 
point distribution. See 3.1.2. 

DNSYH 

Integer which, when $0, indicates an unsymmetric 
wing. 


EDI 

Number of discontinuities 
edge. If there are none, 
1 and a fake kink will be 
if LITYPE-0. 


(kinks) on the leading 
then the number will be 
put in (YDL) . NDL-HL 


122 

Humber of kinks on the trailing edge. This 
number has a minimum value of 1. NDT=NT if 
TETYPE=0. 
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GEOMETRY PROGS A 8 


M 

Reference number for spanwise control points. 
See 3.1.20. In subsequent programs this value is 
referred to by the name MREF. 

JJHAX 

The maximum number of integration stations. 

MFI APS 

The number of flaps. 

(TITLE <26 )1 

Alphanumeric titling information. See 3.1.1 
NTITL 

The ” number cf words in the title. 
NTITL=80/ (integer word length of the computer 
being used) . 


3.2.2 5ECCHD LOGICAL RECORD 


(CHI CP (PP) ) 

Chcrdwise control points. 

( NIBDEX (NF) ) 

Integer array frcm which the spanwise control 
points may be derived using (ETA). 
ETACP (I) =ETA (NTNDEX (I) * (JJHAX+ 1 ) / (NN 4- 1) ) . See 

below for (ETA) . 

(TABLE T (JJHAX ) ) 

Tangents of the sweep angles of the leading edge 
just to the left cf the integration points. This 
array is denoted by LAMLEL in the program, but 
LAHLEL is a aisncmer. The next 3 arrays have 
also been renamed. 

(TAB LES (JJHAX ) ) 

Tangents of the sweep angles of the leading edge 
just to the right of the integration points. 

(TANTE1 (JJHAX) ) 

Tangents of the sweep angles cf the trailing edge 
just to the left cf the integration points. 

(TASTES (JJ HAX) ) 
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Tangents of the sweep angles of the trailing edge 
just to the right of the integration points^ 

(ETA { JJWAXn 

This is a double precision array giving the 
spanwise locations of the integration points. 
The values are nonditnensicnalized by the 
effective semispan and are given by the following 
equation: FT ft (J) - CCS (J*PI/ (JJSAX + 1) ) . 

ISIHETfl (JJH AX) ) 

This is a double precisicn array related to the 
(ETA) array. STHETA (J)=SQPT ( 1 . -ETA (J ) **2) . 

(XSIIIT (JJ MAX) ) 

This array gives the streamwise coordinates of 
the leading edge at the integration stations. 
Values are nondiiensionalized by the effective 
semispan. 

(CCBDir (JJ MAX ) ) 

This array contains the chord lengths at the 
integration stations. Values are 

nondimensionalized by the effective semispan. 

BEAT 

lateral reference length/effective semispan. 

CB A BE R 

Reference chord/effective semispan. 

AR 

Aspect ratio. 

TR 

Sing taper ratio. 

MACH 

Default Mach number (floating point). 

( YEL IBD T) ) 

Spanwise locations of kinks in the leading edge 
nondimensionalized by the effective semispan. If 
no kirks actually exist, then there will be 1 
number here anyway ana its value is 1.E30. When 
LET YPE=0 , this array is the same as (YL) except 
for being nondimensicnal (See 3.1.7). 

JXCTJRDT11 
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Trailing edge counterpart of the above. 
3.2.3 THIRD ICGICAL RECORD (NFLAPSYO) 


HACH 

Mach number again. 

LAME AC 

Sweep angle of the hingeline in radians. 

CBFCES 

Ratio of the reference chord of the flap to the 
reference chord of the wing. 

ftjPAW 

Ratio of the reference area of the flap to the 
reference area of the wing. 

Span wise location of left flap edge. 

ETA 2 

Spanwise location of right flap edge. 


AUZJ21L 

The coordinates of the leading edge at the flap 
side edges. If there are two flaps, these values 
apply for the flap on the right side of the wing. 
The values have teen nondimensicnalized by the 
effective semispac. The first value is for the 
leftmost side edge. 

JCYFJ211 

Chord lengths at the flap side edges. 

XSI 1 

Streamwise coordinate of the left corner of the 
flap. 

XSI 2 

Streamwise coordinate cf the right corner of the 
flap. 

1CHI FPI tJJBAn ) 

This is an array giving the location of the flap 
hingeline at the integration stations. The 
values are the usual CI3X values such that 
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0<CHI<1. 

i£2IPJJ 

Values of the function given by eg. S4e of ref. 1 
at the integration stations. 
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GEOMETRY PROGRAM 


3.3 BOUNDARY CONDITION FILE 


The following data is written on the file referenced 
by the integer W8 in the program. The data is 
written in binary form. The writing of the data on 
the file can be suppressed by inputting NOB as 
non-zero (see 3,1.2). The boundary condition file 
gives the sets of right hand sides for the equation 
solving program. 


3.3.1 FIRST RECORD 


ID 

The identification number. 

ID 3 

ID3=0 for this program. Eor the boundary 
condition program this number is assigned. 

( TIT LE (26) ) 

Alphanumeric titling information. 

UNSYM 

NSYH The number of symmetric cases. NSYM= the 
number of -TRUE, values of (DCS) (See 3.1.3). 

NASIM 

The number of antisymmetric cases. NASYM= the 
number of .TRUE, values of (BCAS) (See 3.1.4). 

JBCS1 

JILASl 

FP 

CWTYPE 

(NF is called MM in ether programs.) 

BMP 

Eor a symmetric wing MMP= (NF*1) /2. 

M MPA 
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For a symmetric wing MMPA=NF/2 . For an 
unsymmetric wing, MMP = MMPA = nf. 

NF . 

.SUIT PE 

JCHICjyPPll 

(NINDEX ( H H P ) ^ 

(E1ACP ( EMF)) 

The control point locations given in 3.1.22. 
Note that this array is single precision and it 
is not necessary in view of having also stored 
(NINDEX (MNP) ) , but it was merely found convenient 
to do it this way. Note also that if the wing is 
symmetric, then cnly the control points actually 
used for the symmetric cases are written. 
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4 METHOD OF SOIUT ICN 


The areas and mean geometric chords were computed first 
by making the substitution y=(b/2)cos€ and using 
SUBROUTINE INTGRT. This subroutine is self -documented. 
The approach outlined in SUBROUTINE INTGRT was highly 
desirable because the integration locations for this 
subroutine coincide with the integration stations which 
must be computed and stored anyway. SUBROUTINE INTGRT 
is called twice for each integration to be performed. 
The first call uses approximately half of the available 
integration stations and the second call uses all of 
the integration stations. The two values are compared 
to give an indication of the error involved. 

A method of weighted quadratics is used for the 
interpolation when the edges are curved. The 
interpolation routine, SUBROUTINE CODIM, was taken from 
ref. 2. 

The methcd of obtaining the flap downwash mode is 
documented in ref. 1 and in SUBROUTINE FLPEWN. 
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5 USE OF TRANSFO RM ATI O N C APA BILITY 

The capability to yaw, shear, or rotate the wing 
applies only to the wing edges as defined fcy the arrays 
(YL) , (XL), (YDL), (YT) , (XT) , and (YCT) . If LETYFE=0, 

then (YDL) is identical to (YL) . Otherwise (YDL) is a 
subset of (YL) derived from (THDDL) (see 3.1.10). 
(YDT) is similarly related to (YT) . After (XL) and 
(XT) are multiplied fcy the scale factors LFSCAL and 
TESCAL the transformation is performed. This is all 
that the transformation capability does. The 
transformation is accomplished upon reaching statement 
1190 of the program and just prior to reading ITFM 16 
of the input. All program actions and input beyond 
this refer to the transformed wing. In particular, if 
the program computes the wing and/or flap reference 
chord, the computed value will be that of the 
I transformed wing. This will generally be different 
from the unt ransf ormed wing if the wing has been 
transformed by yawing. 

When TRNFRWCO (i.e., a cranking transformation), the 
leading and trailing edge definitions must include 
Y=0.0 and Y=0.0 becomes a kink which must be named in 
(INDDL) and/or (INDDT) if either or both is read. 

Ihen the wing is to be transformed by rotating, the 
user must define the leading edge as the portion of 
the periphery which becomes the leading edge, not the 
portion of the periphery . which is the leading edge 
prior to rotation. The same is true of the trailing 
edge. For example, suppose the wing were a square 
defined by the coordinate pairs (0,1), (0,-1), (2,-1), 

and (2,1). Ordinarily the leading edge would be 
thought of as the straight line segment joining (0,1) 
and (0,-1) and the trailing edge as the segment joining 
(2,1) and (2,-1). However, if the user were to want to 
yaw such a wing (TRNANOO) , he would have to tell the 
program that the wing is unsymmetric (UNSYMfO) and that, 
the leading edge is comprised of the two line segments 
joining (2,1) with (0,1) and (0,1) with (0,-1). 
Similarly, he would have to tell the program the the 
trailing edge is comprised of the two line segments 
joining (0,-1) with (2,-1) and (2,-1) with (2,1). 

The transformation capability is useful for inputting 
simple, swept, tapered wings when the leading edge 
sweep angle is known. For this situation make TRNFRIKO 
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(see 3.1.2); KL=3 and TBHAHG - leading edge sweep angle 
(3.1.6); TL = semispan (3.1.7); card 8 = blank (3.1.8); 
NT=3 (3.1.11); YT { 1 ) = semispan (3.1.12); and XT ( 1 ) = 
tip chord and XT (2) - root chord (3.1.13). 

The transformation capability is also useful for delta 
wings. For this situation snake TRNFFM<0; NL=3 and 
TBNARG = leading edge sweep angle; YL(1)= semispan; 
card 8 - blank; NT= 1 ; card 10 = blank; and XT(1) = root 
chord. 



GT CHET BY FFOGBAH 


6 SAMPLE CAS ES 


Some sample cases are given in appendix A. Note that 
blank fields are equivalent to zerc for numerical data 
and to .FALSE, for logical data. Fig. 6 and fig. 7 are 
plots of the two sample cases. These plots were made 
with a planform plotting program which uses the 
geometry file created by the geometry program. 
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7 CONT ROL CARPS 

Job centrcl cards needed to handle the files are given 
below. 


7.1 


Ames* IB M T5S system 

For the first time that the program is to be used 
the following TSS commands must be given. These 
need to be given once and only once for each user 
ID. 

SHARE MEDAN, FS ARTB , I Nl DFILF 
CCS MEEAN , IDFILE 
DELETE MEDAN 

SHARE MEDAN, FS ARTE, ISP BCG. V 1 

After the above commands have been issued by the 
user, the program is invoked by the following, 

JBLB MEDAN 
CALL GECHS 
(input data) 

For the Ames* TSS computer system no control 
cards are needed for the geometry and boundary 
condition files. The DDEF cards are 

automatically issued by the GEMFIL and BCFIL 
subroutines for the . geometry file and the 
boundary condition file, respectively. The DDEF 
card for the identification number file is issued 
directly from the main program using the Ames' 
library programs CVRT and OBEY. The latter is 
the one which actually issues the control card tc 
the operating system. The DSNAMES (data file 
names) assigned tc the geometry file and boundary 
condition file would be GEOM.X23 and BC.X023.XO00 
if ID1 were equal to 23. The DSNAME of the 
identification number file is always IDFILE. 

Because of this automatic file defining feature, 
the following messages may be received in the 
output from the program. These messages are 
normal and do not indicate any errors. 

CANCELLED: DDNAME ET07F001 UNKNOWN 
CANCELLED: DDNAME FTOBElOl UNKNOWN 


26 



G1CHETBY PBOGBAB 


CANCELLED : DDNAHE F109F301 UNKBCWN 

The program reads from unit 5 and writes on unit 
6, so BEEF commands such as the following may be 
required if the program is run conversationally, 

EDEF FTQ5F001,, input 
BDEF FTC6F001, , output 
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JHPEKDIX A 





FIRST SAMPLE CASE 





c o 1 u rou 



1 

2 

3 4 5 6 7 R 

HEM 

123456789012345678901234567890 12 34567890 12345678 90 1 2 34567 R 90 1 2 34567 8 90 1 234567890 

1 


SAMPLE SING 

1 

2 

0 

1 0 

1 0 

0 0 2 0 0 

3 

FFEFFFFFFF 



4 

FF FFFFFFFF 



e 

0.0 




6 

7 

1 

0.0 


7 

1 . 

0.7 

0.5 

0.0 

8 

0.5 

0.5 

0.5 

0 . 0 

9 

2 




10 

3 




11 

5 

1 



£ 12 

1.0 

C.7 

0 . 0 


13 

1.5 

1.5 

2.0 


14 

0 




16 

0 . 

0 . 

0 . 


19 

71 




20 

5 

17 





SECOND SAMPLE CASE 







column 





1 

2 

3 

4 5 

6 

7 8 

ITEM 

12.34567890 12345678901234567890 12345678 9012 345673901 2 34567890 1 23 45678 9012 34567390 

1 


SAMPLE SING 

2 




2 


1 


1 1 

2 



3 

FFFFFFFFFF 






4 

FFFFFFFFFF 






5 

C.O 







6 

9 

1 

.5 

20.0 




7 

1 . 

1 . 

1 . 

C. 7 

0.5 0.0 

-0.5 

-0.7 

1 

- 1.0 







9 

1.5 

1.0 

0.5 

0.5 

0.5 0.0 

0.5 

0.5 

9 

0.5 







9 

3 







10 

3 

5 7 






11 

7 

1 






12 

1.0 

e. 7 

0 . 0 

-0.7 

- 1.0 - 1.0 

- 1.0 


13 

1. 5 

1.5 

2.0 

1.5 

1.5 1.0 

0.5 



15 5 

16 

19 10000 

20 5 
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■ Hinge line extension at 
constant % chord 


Figure 1.- Definition of coordinate system. 
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Original Wing 


TRNFRM=0 



TRNFRM <0 


TRNFRMX) 


Figure 5.- Illustration of the transformation capability for TRNANG = 45 
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C0NTR0L P0INT L0CRTI0N 

WINE 1RTR 
R5PEET RHTIH 
LQNG./LRT. REF. LENGTH 
TRPER RRTI0 



ii ii n 

1.53527 

1.39862 

.£□□□□ 


5RMPLE WINE 

1 



Figure 

6.- The planform and control 

points 

of the 

first 

sample case. 



2D □ 


V 



WINE HATH 
fiSPECT RHTIB 
LQNE./LAT. REF. LENETH 
TAPER RATID 


+ CONTROL POINT L0CRTI0N 

= 1.S3D4-7 

= 1.3S6D1 

= .□□□□□ 


5RMPLE WINE 2 

Figure 7.- The planform and control points of the second sample case. 
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FLAPPED WING GEOMETRY PROGRAM 


TITLE OF CONFIGURATION = SAMPLE WING I 

ASSIGNED IDENTIFICATION NUMBER = 19 


C W TYPE * 1 
SWTYPE = 0 
I P INFO = 1 
CPRESS = 0 
UNSYM = 0 
TRNFRM = 0 
OUT = 2 
NOT = 0 
N08 » 0 



INTEGRATIONS FOR THE WING AREA 

VALUE OF THE INTEGRAL IS 2.6042147 THE NUMBER OF INTEGRATION POINTS USED = 35 

NONE OF THE X-VALUES WERE CONSIDERED DISCONTINUITIES 
THE X VECTOR IS 

0.500000 -0 .500000 

VALUE OF THE INTEGRAL IS 2.6053963 THE NUMBER OF INTEGRATION POINTS USED = 71 


INTEGRATIONS FOR THE WING CBAR 

VALUE OF THE INTEGRAL IS 3.6425514 THE NUMB FR OF INTEGRATION POINTS USED = 35 

NONE OF THE X-VALUES WERE CONSIDERED DISCONTINUITIES 
THE X VECTOR IS 

0.500000 -0.500000 

VALUE OF THE INTEGRAL IS 3.6439657 THE NUMBER OF INTEGRATION POINTS USED = 71 

u> 



GEOMETRICAL DATA FOR THE FLAPPED WING 


WING DATA 


WING AREA = 2.60540 
BREF, LAT. REF. LENGTH = 1.00000 
B2t EFFECTIVE SEMI-SPAN = 1.00000 
BREF/B2 = 1.00000 
CBAR » LONG. REF. LENGTH = 1.39862 
ROOT CHORD = 2.00000 


ASPECT RATIO 
TAPER RATIO 
BREF/CBAR 
CBAR/ROOT CHORD 
MACH NUMBER 


1.53527 
0. 50000 
0.71499 
0.69931 
0.00000 


THE NUMBERS IN PARENTHESES ARE ERROR ESTIMATES 
(SEE PROGRAM DOCUMENTATION) 


( 0.00118) 

( 0.00054) 

< 0.00070) 

( 0*00028 ) 

( 0.00027) 



LEADING EDGE CO-ORDINATES 


YL 

1 .OOOOO 
0.70000 
0.50000 
- 0.00000 
-0.50000 
-0.70000 
-1.00000 


TRAIL IN G EDGE 
YT 

1.00000 

0.70000 

0.00000 

-0.70000 

- 1.00000 


XL 

0.50000 

0.50000 

0.50000 

0.00000 

0.50000 

0.50000 

0.50000 


CO-ORDINATES 

XT 

1.50000 

1.50000 

2.00000 

1.50000 

1.50000 


(SCALE FACTOR = 1.000000) 


(SCALE FACTOR = 1.000000) 



THE WING HAS KINKS IN THE LEADING EDGE AT THE FOLLOWING ETA LOCATIONS 


0.50000 

- 0.50000 



GEOMETRICAL DATA AT THE INTEGRATION STATIONS 


ETA 

X ( LE>/82 

C / B2 

T A ML EL 

TANLER 

TANTEL 

TANTER 

0.99905 

0.50000 

1.00000 

0.00000 

- 0 .00006 

0.00000 

- 0.00100 

0.99619 

0.50000 

1 .00000 

0.00000 

- o.ooooo 

0.00000 

- 0.00000 

0.99144 

0.50000 

1.00000 

0.00000 

- 0.00000 

0.00000 

- 0 .00000 

0.9 84 81 

0.50000 

1.00000 

“ 0.00001 

- 0.00000 

0.00000 

- 0.00000 

0.97630 

0.50000 

1.00000 

0.00000 

0.00001 

0.00000 

- 0.00000 

0.96593 

0.50000 

1.00000 

0.00000 

- 0.00000 

0.00000 

- 0.00000 

0.95372 

0.50000 

1.00000 

0.00000 

- 0.00000 

0.00000 

- 0.00000 

0.939 69 

0.50000 

1.00000 

“ 0.00000 

- 0.00000 

0.00000 

- 0.00000 

0.92388 

0.50000 

1 . 00000 

0.00000 

0.00000 

0.00000 

- 0.00000 

0.90631 

0.50000 

1.00000 

0.00000 

- 0.00000 

0.00000 

- 0.00000 

0.88701 

0.50000 

1.00000 

0.00000 

- 0.00000 

0.00000 

- 0.00000 

0.8 6603 

0.50000 

1.00000 

- 0.00000 

- 0.00000 

- 0.00004 

- 0.00000 

0. 84339 

0.50000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00004 

0.81915 

0.50000 

1.00000 

0.00000 

- 0.00000 

0.00003 

- 0.00000 

0.79335 

0.50000 

1.00000 

- 0.00000 

- 0.00000 

0.00000 

- 0.00003 

0.76604 

0.50000 

1.00000 

0.00000 

0.00000 

0.00000 

- 0.00000 

0.73728 

0.50000 

1.00000 

0.00000 

- 0.00000 

- 0.00003 

- 0.00000 

0.70711 

0.50000 

1.00000 

0.00000 

- 0.00000 

- 0.20180 

0.00003 

0. 67559 

0.50000 

1.00 636 

- 0. 00000 

- 0.00000 

- 0.36513 

0.20180 

0.64279 

0.50000 

1.01834 

0.00000 

0.00000 

- 0.47922 

0.36513 

0 . 60876 

0.50000 

1.03464 

- 0.00000 

“ 0 .00000 

- 0.58765 

0.47922 

0.57358 

0.50000 

1.05532 

- 0.00002 

0.00000 

- 0.68845 

0.58765 

0.53730 

0.50000 

1.08029 

0.00002 

0.00002 

- 0.77967 

0 . 6R845 

0.50000 

0.50000 

1.10937 

0.42880 

- 0 . 0000 ? 

- 0.85912 

0.77967 

0.46175 

0.4 83 60 

1.15864 

0.69587 

- 0.42880 

- 0.92452 

0.85912 

0.422 62 

0.45637 

1.22204 

0 . 9315 ? 

-0.69587 

- 0.97328 

0.92452 

0.38268 

0.41917 

1.29811 

1.12920 

- 0.93152 

- 1.00278 

0.97328 

0.34202 

0.37325 

1.38480 

1.28162 

- 1.12920 

- 1.01004 

1.00278 

0.30071 

0.32030 

1.47948 

1.38072 

- 1.28162 

- 0.99195 

1.01004 

0.25882 

0.26247 

1. 57886 

1.41752 

-1 .38072 

- 0.94514 

0.99195 

0.21644 

0.20239 

1 . 67899 

1.38202 

- 1 .41752 

- 0.86601 

0.94514 

0.17365 

< 5 . 14326 

1.77519 

1.26304 

- 1.38202 

- 0.75061 

0.86601 

0.13053 

0.08879 

1.86202 

1.04805 

- 1.26304 

- 0.59474 

0.75061 

0.08716 

0.04334 

1.93326 

0.72297 

- 1.04805 

- 0.39392 

0. 59474 

0.043 62 

0.01186 

1.98189 

0.27190 

- 0.72297 

- 0.14332 

0.39 392 

- 0.00000 

0.00000 

2.00000 

- 0.27190 

- 0.27190 

0.14332 

0. 14332 



CHORDWISE CONTROL POINT LOCATIONS REFERENCED TO THE LOCAL CHORD 


0*07937 

0.29229 

0*57116 

0.82743 

0.97975 


SPANWISE CONTROL POINT LOCATIONS REFERENCED TO B2 


0.98481 
0.93969 
0.86603 
0.7 6604 
0. 64279 
0.50000 
0.34202 
0 . 1 73 65 
- 0.00000 
-0.173 65 
-0.34202 
-0.50000 
-0. 64279 
-0.7 6604 
-0.86603 
-0.93969 
-0.98481 


= 71 


JJMAX 



THE GEOMETRY INFORMATION HAS BEEN STORED ON UNIT W7 

V 


•P* 

Ln 



NO BOUNDARY CONDITIONS WERE STORED ON UNIT W8 



FLAPPER WING GEOMETRY PROGRAM 


TITLE OF CONFIGURATION = 

ASSIGNED IDENTIFICATION NUMBER = 21 


CWTYPE 


1 

SW TYPE 


0 

I P IN FO 

.... 

0 

C PR ESS 


o 

UNSYM 

= ■ ' 

1 

TRNFRM 


I 

OUT ; 

‘ 

2 

N07 

S 

0 

NOB 


0 


SAMPLE WING 2 




INTEGRATIONS FOR THE WING AREA 


THE UNCORRECTED VALUE OF THE INTEGRAL = 
THE CORRECTED VALUE OF THE INTEGRAL = 
THE NUMBER OF INTEGRATION POINTS WAS = 

ORIGINAL DISCONTINUITIES 

0.692069 0.26905 2 


REDUCED DISCONTINUITIES 

0.692069 0. 269052 

THE UNCORRECTED VALUE OF THE INTEGRAL = 
THE CORRECTED VALUE OF THE INTEGRAL = 
THE NUMBER OF INTEGRATION POINTS WAS = 

ORIGINAL DISCONTINUITIES 

0.692069 0.269052 


REDUCED DISCONTINUITIES 
0.692069 

S 


0.269052 


INTEGRATIONS FOR THE WING CBAR 

THE UNCORRECTED VALUE OF THE INTEGRAL = 
THE CORRECTED VALUE OF THE INTEGRAL = 
THE NUMBER OF INTEGRATION POINTS WAS = 

ORIGINAL DISCONTINUITIES 

0.692069 0.269052 


REDUCED DISCONTINUITIES 

0.692069 0.269052 

THE UNCORRECTED VALUE OF THE INTEGRAL = 

THE CORRECTED VALUE OF THE INTEGRAL = 

THE NUMBER OF INTEGRATION POINTS WAS = 

ORIGINAL DISCONTINUITIES 

0.692069 0.269052 


REDUCED DISCONTINUITIES 


2. 1186IB0 
2. 1185694 
143 


-0.5769R3 -0.692069 


-0.576983 -0.692069 

2.1185522 
2. 1185484 
287 


-0.576983 -0.692069 


-0.576983 -0.692069 


2.6438274 

2.6436501 

143 


-0.576983 -0.692069 


-0.576983 -0.692069 

2.6436939 

2.6436787 

287 


-0.576983 


-0.692069 



690269 * 0 - 


£ 8692 . 9 * 0 - 


290692*0 


690269 *0 



GEOMETRICAL DATA FOR THE FLAPPED WING 


WING DATA 


WING AREA = 2.61357 
BREF, LAT. REF. LENGTH = 1.00000 
B2» EFFECTIVE SEMI-SPAN = 1.11070 
RREF/B2 = 0.90033 
CBARt LONG. REF. LENGTH = 1.38601 
ROOT CHORD = 1.64005 


ASPECT RATIO 
TAPER RATIO 
BREF/CBAR 
CBAR/ROOT CHORD 
MACH NUMBER 


1.53047 

0.00000 

0.72149 

0.84510 

0.00000 


THE NUMBERS IN PARENTHESES ARE ERROR ESTIMATES 
(SEE PROGRAM DOCUMENTATION) 


( 0.00003) 

( 0 . 00001 ) 

( 0 . 0000 ?) 

( 0 . 00001 ) 

( 0 . 00001 ) 



LEADING EDGE 
YL 

I. 45272 
1.28171 
1.11070 
0.82879 
0. 64086 
0.00000 
-0.29884 
-0.48677 
-0.76868 


TRAILING EDGE 
YT 

1.45272 
1.17081 
0. 6R404 
•-0. 14475 
-0.42666 
-0.59767 
-0.7 68 68 


CO-ORDINATES 

XL 

1.06752 
0.59767 
0. 12783 
0.23043 
0.29884 
0.00000 
0. 64086 
0.70926 
0.81187 


CO-ORDINATES 

XT 

1.06752 
1.17012 
1.87938 
I .64895 
1.75156 
- 1.28171 
0.81187 


(SCALE FACTOR = 1.000000) 


(SCALE FACTOR = 1.000000) 



THE WING HAS KINKS IN THE LEADING EDGE AT THE FOLLOWING ETA LOCATIONS 


1.11070 

0.64086 

-0.29884 


THE WING HAS KINKS IN THE TRAILING EDGE AT THE FOLLOWING ETA LOCATIONS 


”0.42666 

THE ABOVE VALUES WERE DETERMINED BY APPLYING 
A 20.000 DEGREE ROTATION TRANSFORMATION TO THE INPUT 



CHORDWISE CONTROL POINT LOCATIONS REFERENCED TO THE LOCAL CHORD 


0.07937 

0.29229 

0.57116 

0.82743 

0.97975 


SPANWISE CONTROL POINT LOCATIONS REFERENCED TO R2 


0.98481 
0.93969 
0.86603 
0.7 6604 
0. 64279 
0.50000 
0.34202 
0.17365 
- 0.00000 
-0.17365 
-0.34202 
-0.50000 
-0. 64279 
-0.76604 
-0.86603 
-0.93969 
-0.98481 


JJMAX = 287 



THE GEOMETRY INFORMATION HAS BEEN STORED ON UNIT W7 



NO BOUNDARY CONDITIONS WERE STORED ON UNIT W8 



■n n n n rv o n r* r* o nr* nr* nr* on r* n >r* nn-nnno on 


c 

c LIFTING SURFACE GEOMETRY PROGRAM 
C 

C THE WING OUTLINE CAN BE 

C DEFINED BY UP TO 100 BREAKPOINTS ON THE LEADING EDGE AND UP TO 
C 100 BREAKPOINTS ON THE trailing EDGE 

SUBROUTINES USED 

0 A FUNCTION SUBPROGRAM USED TO CALCULATE A QUANTITY ASSO- 

CIATED WITH THE FLAP PRESSURE MODE. 

C2I A FUNCTION SUBPROGRAM USED TO CALCULATE THE FUNCTION 

GIVEN BY EG, AUE OF NASA TN D-7251 (1973), 

OBEY AN AMES' LIBRARY PROGRAM TO PASS COMMANDS TO THE TSS 
OPERATING SYSTEM", OrEY IS USEO TO release AND ODEF 
UNIT 9, WHICH IS THE IDENTIFICATION NUMBER FILE, 

CODIm A SUBROUTINE SUBPROGRAM USED FOR INTERPOLATION, SEE 
PROGRAM DOCUMENTATION. 

FLPDwn a SUBROUTINE SUBPROGRAM USED TO DETERMINE THE CONTROL 
SURFACE DOWNWASH MODE, SEE PROGRAM DOCUMENTATION, 

STRAIT A SUBROUTINE SUBPROGRAM USED FOR STRAIGHT LI N E INTERPOLA- 
TION. 

INTf,RT A SUBROUTINE SUBPROGRAM used FOR NUMERICALLY INTE- 
GRATING TO FIND THE WING AND CONTROL SURFACE AREAS 
AND mean geometric CHORDS, 

gemfil a subroutine subprogram which issues a release and doef 

command to THE AMES' TSS OPERATING SYSTEM FOR THf 
GEOMETRY FILE, GEMFIL WOULD NOT BE NEEDFD ON ANOTHER 
COMPUTER SYSTEM, 

BCFIL A SUBROUTINE SUBPROGRAM which ISSUES A release AND ddef 
C0MM 4 nd FOR THE BOUNDARY CONDITION FILE, BCFIL WOULD 
WOT flE NEED ON ANOTHER system. 


c SPECIFICATION STATEMENTS 

c 

COMMON f IRSTfXiCH^SQX.DcP 
COMMON ETALI 



n n ri 


COMMON ETATI 
COMMON XDIS 
COMMON XL 
COMMON xi£ 
COMMON x un- 
common XT 
common XTE 
COMMON XTP 
COMMON VOIS 
COMMON YU 
COMMON YUP 

common yt 

COMMON 20 1 S 
COMMON INOOU 
common indot 


double precision 01 

DOUBLE PRECISION eta 
double PRECISION ETAIP 
DOUBLE PRECISION ETaZ 
DOUBLE PRECISION DTHETa 
DOUBLE PRECISION PIO 
DOUBLE PRECISION STHETa 
double PRECISION TMETaD 
INTEGER CPRESS 
INTEGER C*TYPE 
INTEGER OUT, n'JT? 

integer p 

INTEGER TRNFRM 

integer pp 

INTEGER PPl 

integer rs 
integer 

INTEGER SNTYPE 
INTEGER tetype 
INTEGER TITLF 
INTEGER unsym 



INTEGER 
INTEGER to 7 
INTEGER ^8 


LOGICAL BCiS 
LOGICAL 0CS 
LOGICAL LOCI 
LOGICAL PAIR 
LOGICAL TPU M 
C 

real lamdac 
real LAMDCD 
REAL LAMLEL 
REAL LAMLER 
REAL L ANT EL 
REAL LANTER 
REAL LESCAl, 

real linit 

REAL naCH 
C 

DIMENSION ALFAB (15, U7) 
DIMENSION PCS ( 10) 
or DIMENSION cHI C P C 1 S > 

DIMENSION COROFC 383) 
DIMENSION cTF(2) 
DIMENSION ETa( 383) 
DIMENSION ETALI(IOO) 
DIMENSION LAMLER ( 383) 
DIMENSION , N INDEX (47 ) 
DIMENSION TITLE(2fe) 
DIMENSION XLflOO) 
DIMENSION XLF C 2) 
DIMENSION XTF( 393) 
DIMENSION YOL(IOU) 

dimension YLF (?) 
DIMENSION XD| (104) 
DIMENSION INDDT(lOa) 


,ALPHA(4,15,47),9CA$(10) 
,CESG( 383) f CHI (50) 
,CHIFPI( 383) 

, CORD I P ( 383) , C ORDSQ ( 383) 

,C2IP( 383) , OCR ( 50 ) 

i ET aF (2, 8) ,ET4IF(2) 

,ETATI(100) » LAMLEL ( 383) 

,LAMTEL( 383) ,LAMTEH( 383) 
, SOX f 5 0 ) ,STHETA( 383) 

, X (50 3 ,XDTS(200) 

, XLE ( 383) 

# X S I L I p ( 383) #XT(iOO) 

, X TF ( 2 ) ,YDIS(200) 

,YOT(10a) , YL ( 100) 

,Yt(100) ,ZDIS(200) 

#XDT<104) ,INODlMO«) 

,ETACP(47) 


c 

c 
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c 

EQUIVALENCE fLAMLER#CORDSQ) 

EQUIVALENCE (LAHTELf COffOF) 

EQUIVALENCE (LAMTf.R,CFSQ) 

C 

C 

C LINEAR INTERPOLATION FUNCTION FOR X»F(Y) 

e 

LINIT(V,Xi,X? r Yi f Y 2 )»Xi+(V-Yn/ ( Y 2 *Yl)*(X 2 .Xi) 


DU1.D0 
ETAZsO'.DO 
JJMXMXs 3 B 3 

PI *3, 141593 

PtD»3, 141592853589793 00 

R2Ds57. 29578 

BETAPsO, 

beta*o; 

TANSHR so, 

R5=5 

W9*9 

Wfc *6 

W7»7 

we«8 
YAvGsO , 

VDLtl)=l,E30 
YDT ( 1 ) * 1 , E30 

, I R 0 R P L SHOULD pE SET TO THE INTEGER RQRD LENGTH ( I N CHARACTERS) 

,0F TwE PARTICULAR COMPUTER RE I Nr, USED TF TmE PLOT PROGRAMS ARt TO RE 

.used, 

JwOROLm 


TITLE ■ HEADING FOR PRINTED OUTPUT 

ID s IDENTIFICATION NUMBER 

CRTYPE • TYPE OF CHORD«ISE CONTROL POINTS 
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O' 

o 


?5 

C 

c 


c 

c 

c 

c 

c 


SWTYPE 

IPINFO 

CPRESS 

UNSYm 

TRNFRM 

OUT 

NOT 

NO? 


TYPE OP SPANWISE CONTROL POINTS 

PRINT INFORMATION at INTEGRATION STATIONS 

COMPUTE delta CP(FLAP) at nc chordrise control STATIONS 

UNSYmETRICAu wing Pi aNFORM lut . BIU . Tu _ UTwr , 

INTEGER FLAG FOR YaWING, CRANKING, OR SHEARING THE WING. 

OUTPUT LEVEL FROM SUBROUTINE INT(;RT 
FLAG TO PREVENT WRITING ON UNIT W? 

FLAG TO PREVENT WRITING ON UNIT Wg 


NTITL * SO/IWOROL 
READ (RSfSOOt) (TITLE(N) . N * 
RF AD (95, 5003) ID, CWTYPe, 
l OUT, NO?, NOS 
0UT2 a MINO<OUT#2) 
IFMD.NI.O) GO TO 25 


1# NTITL) 
SWJYPf, ipinfo, 


CPRfSS, UNSYM* TRNpRM, 


FILE DEFINITION for aMES» TSS SYSTEM, 
call OBEV(U,UHRELEaSE FT09F001 ) 

CALL OBEY (22,2 ttHODEF FTQ9F001 , , IOFI^E 
REWIND 9 

REAP ( W9) ID, I D i , IDS. ID3 
lORlDtl 
REWIND 9 


WRITE ( W 9 ) ID, 
ENp PILE WR 
continue 


IDI, ID2, ID3 


WRITE (WS, >01 9) 

WRIT F(Wfe, 6006) (TITLE (N) , N* j , NT iTL) 

WR I f E ( W6 j 6021 ) ^wTYPE, SwTYPE, IPINFO, CPRESS, UN5YM , TRNFRM, OUT,. 
l BC$ 7 ' N ' 08 * LOGICAL Array SPECIFING the SYMETpIC BOUNOARY CONDITIONS 

READ (R5,500«) (BC8(1 ) ,1*1,10) 

BCAS a LOGICAL ARRAY SPECIFIING THE ANT I -SYMMETR IC BOUNDARY CONdITIOK 
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RfcAD (R5,S00«) (BCASrX)'lBl, 10 ) 

NF|. APS ® 0 

IP (8CS(«},0R,8CA5(3)) NFLAPS =-l 
MACH s MACH NUMBER 

Re ao (R5» S 0021 mach 
109 CONTINUE 

N i * number of co-ordinate pairs for leading edge definition 

LETYPE s TYPF OF INTERPOLATION USED ON LEADING EDGE 
LESCAL s SCALE FACTOR XL ( ! ) »XL ( I ) *LESC AL 

XcLt * INTERPOLATION CONTROL CONSTANT FOR end INTERVALS 

AND LETvPE ,NE, 0 

TRNANG 8 SIDESLIP, CRANK, OR shear angle in degrees 

read CR5, 5005) NL , LETYPE, LESCAL, XCLE, TRNANG 
IF(NL,LT,2) NL*2 
IF(LE5CAL»EG.0,0) LESCALBl.O 
IF(UNSYM # EQ,0) NLPs fNL + n/2 
TF ( UNS VM f Nf , 0 ) NLP.NL 

Yid) B SPAN wJSf CO-ORDINATES OF LEADING EDGE 
XLM) s CHORDHISE CO-ORDINATES OF LEADING EDGE 

READ (R5,5002) ( YL f I ) , I « ) , NLP ) 

READ (R5# 500?) (XL(I),I«1,NLP) 

IF ( UN$ Ym , NE , 0 ) GO To 1102 
C IF (UNSYM.EQ.O) CONTINUE 

DO 1101 I c 1, NLP 
NLL * NL ♦ 1 » ! 

YL (NLL) » -YL ( ! ) 

XL f NLL ) s XLU) 

1101 CONTINUE 

1102 CONTINUE 

BREED > (YL ( 1 )*YL(NL) )/2, 

00 110 I s 1, NL 

*L< I)»XL(I)*LE8CAL 

no CONTINUE 
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NDL a NUMBEH of DISCONTINUITIES on LEADING edge 

VOLCI) * V LOCATION. OF DISCONTINUITIES 

INDDL « INTEGERS FOR DETERMINING YdL when LETYPE ,NE, o 

IF (LETYPE.NE, 0) GO TO 1104 
NDL * Nl 

DO 1103 I a 1, NDL 
YDL C I > 5 y U C I > 

XOL ( I ) * XL(I) 

1103 CONTINUE 

GO to mo 

1104 continue 

HEAD (R5,5003) NDL 

IF (NpL ,£Q. 0) GO TO 1110 

IF (UNSYm* N£, 0 ) GO To 1107 

IF (UNSYH.EO.03 CONTINUE 

NOLP * (NDL* 1 3/2 

NOLI « NOL *1 ■ „ 

READ (R5, 5003) ( ZNDDL, f 1 5 # I*l» ^PLP) 

DO 11061 a 1, NOtP 
NDLE * NDL 1 ~ I 

INDDL(NDLE) a NL * 1 - INDDL(I) 

1106 CONTINUE 
GO TO 1100 

1107 CONTINUE 

NDL p * NDL 

READ (R5, 5003) (INDDL(I), 1*1# NDL p ) 

1106 CONTINUE 

DO 1109 I a 1, NDL 
IND a INDOL 
XDL(I) a XL(IND) 

YDL(I) a YL(IND) 

1 1 OR CONTINUE 
1110 CONTINUE 

NT a NUMBER OF CO»ORDINaTE PAIRS FOR TRAILING EDGE QEFINITIO 

TETYPE a TYRE OF INTERPOLATION USED ON TRAILING EDGE 

TESCAL a SCALE FACTOR XT(I)aXT(I)*TESCAL 
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XCTE * interpolation CONTROL constant for END INTERVALS 
AND TETYPF ,NE. 0, 

READ (R5,5005) NT , TET YPE , TESC aL , XCTE 
IF fNT.LT, 2) NTb 2 
IF(UNSYM,EQ,0) NTPe (NT<M)/2 
IF fUNSYN.NE.O) NTPaNT 
IF fTESCAL.EO.O, O') TESCALeLESCAL 

YT m a SPAN* I SE CO-ORDINATES OF TRAILING EDGE 
X T ( I ) a CHORDWISE CO-ORDINaTFS OF TRAILING EDGE 

READ (R5,b002) C YT f I ) , 1*1, NTP) 

READ <P5/500*) (XT(I),I*WNTP) 

If (Unsym.ne.O) go to i 1 \i 

IF (UNSYM.EN.O) CONTINUE 
NT A s nTP . 1 

nTB * wT * \ 

DO ilii I « 1, NTA 
NTC a NTB . I 

YT(NTC) 8 wYT(T) 

XT(nT C ) s XT(I) 

1111 CONTINUE 

1112 CONTINUE 

DO ill I s 1 | NL 
XTtI)8XTri)*TtSCAL 
til CONTINUE 

NDT s NUMBER OF DISCONTINUITIES On TRAILING EDGE 

YDT(I) e Y LOCATION OF TRAILING EDGE DISCONTINUITIES 

IF ( TFT YPE , Nf t o ) GO TO lll« 

NDT s mT 

DO 1113 I a 1 # NDT 
YDT (I } c VT(I) 
x o t ( i ) s xtci) 

1113 CONTINUE 
GO TO 1120 

1114 CONTINUE 
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RE AD (R5# 5003) NOT 
IFfMDT t EG, 0) GO TO 1120 
IF CUNSYm.ne.O) GO TO 1117 
C If (UNSYM.EC.O) CONTINUE 

NpTP s CNOT+i )/? 

NOT | s NOT ♦ 1 

RE AH ( R5f 5003) fINDDTCl)# I«lf N&TP) 

DO 1116 I * l, NDTP 

NOTE ■ NOT 1 - I 

INODT(NDTE) s NT ♦ 1 - INODTC!) 

1116 CONTINUE 
GO TO 1118 

1117 CONTINUE 

NDTP s NOT 

READ <R5> 5003) ( INDD-T C I ) # U t r NDTP) 

1118 CONTINUE 

DO 1119 I * 1, NOT 
INO a INDDT(I) 

XoTCl) * XT (INO) 

. vDT(I) « VT(IND) 

1119 CONTINUE 

1120 CONTINUE 
C 

IF (TRNANG .EG, 0,) GO TO 1190 

.....the wing will be transformed bv cranking, shearing, or rotating. 

9 w w m w v 

IF (TRNERM) 1121, 11 2a, 1 124 

1121 CONTINUE 

C , , , , .THE CRANKING TRANSFORMATION WILL BE APPLIED TO THE X-VALUES, 
TAnShR s TAN(TRNANr/R2D) 

00 1122 I * 1, NL 

XL { I ) * XL ( I ) ♦ ABS(YL(I))*TANSHR 
XOL(I) 8 XoL(I)*AeRCYDLCI) )*TANSHR 

1122 CONTINUE 

DO 1123 I * 1, NT 

XT f I ) e XT ( I ) ♦ A05 ( YT (I) ) * TANSHR 
XDT(I)8XnT(I)TABS(YDTCl))*TANSHW 

1123 CONTINUE 
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GO TO 1190 
112« CONTINUE 

IF (UNSYM ,NF f 0) GO TO 1 125 

it is assumed that the *ing will become unsymmetric if 

• • t • » I T IS SHEARED OR ROTATED'. 

UNSYM » fl88 
NIP s NL 
NTP s NT 
NDLP # NOL 
NpTP s NDT 
1125 CONTINUE 

IF (TRNFRM) 1J50# 1130, 1160 

1130 CONTINUE 

the shearing transformation is applied to the various x-values, 

TANSHR b TAN(TRNANG/R2D) 

00 1135 I » 1, NL 
YK I ) s X L ( I ) - YLTI)*TANShR 
XOL(I)»XDL(I)"YDL(I)*Tanshr 
1135 CONTINUE 

00 1H0 I s 1, mT 
XT (I ) b X?(I) - YTf I)*TANSHR 
XDT (I)bXO'T fI)-YDT (I)aTANSMR 
1UO CONTINUE 
GO TO 1190 
116f) CONTINUE 
C 

c..,,,the rotation by angle trnang is applied to (X,Y) pairs, 

BETA b TRNANG 
BETAR b BETA/R20 
CR b COS(BETAR) 

Sr * SIN(rETAR) 

DO 1165 I b t, NL 
XTEHP e XLfI) 

XL f I ) a C6*XTEMP • SR*¥L(I> 

YL f I ) * SB*XTEMP ♦ C8*YL(D 
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1165 CONTINUE 

DO 1170 I a l, NT 
XTEMP * XT ( I ) 

X T ( I ) a CB*XTEMP - Sa*YT(I) 

YT(I) * SB* xTfcMP + XB*YTCI) 

1170 CONTINUE 

IF (NDL) 1177, H77, 1174 

1174 DO 1175 1 = 1 , NOt 
XTEMP a XOLCI) 

XDLU ) = CB*XTEMP - SB* YDL C I ) 
YDUI) a SB* XTEMP ♦ CB*YDl(!) 

1175 CONTINUE 

1177 IF (NOT) 1 19 q, 1 190 , 1 179 

1179 00 1180 I « 1 , NOT 
XTEMP a XOT(I) 

XDT(T) a CB*XTEMP • Sa*YOT(I) 

' ' YC T ( j ) a XTEMP*SB ♦ CB*YOT(j) 

1180 CONTINUE 
U9o CONTINUE 

Y A V G a (YL( 1 ) ♦ YUNt ))/?,0 

IF(UNSYM, FQ, 0 ) GO TO 113 
8*Y L (n-YL(NL ) 

B ? a B / 2 , 0 
GO TO 114 

113 CONTINUE 
B 2 « Y l ( 1 ) 

a 8.? , 0 * ft 2 

114 continue 

DO H5 Isi.NL 

pTALKI) S (Yim - YAVGl /82 

115 CONTINUE 

DCf 1 1 6 lal , NT 

rTATKI) s (YT(I) - V*Vfi)/B2 

116 CONTINUE 


ARE AN a AREA OF TH£ WING 

C6ARW a mfaN AERODYNAMIC CHORD (REFERENCE CHORD) 
pjREF a STREAMWlSf RING REFERENCE LENGTH 
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AREAF a AREA OF THE FLAP 

C0ARF a CHORD REFERENCE LENGTH OF THE FLAP 

READ (R5,5002) AREA*, CBaRW, BREF, AREAF, CBARF 
IF (BREF ,EQ, 0.) BREF a BREFD 
IF (NFLaPSjEG.O) GO TO 112 

Y F 1 , Y F 2 a SPANWISE POSITIONS OF THE FLAP EDGES (YFl.LT.YF2) 

CHI 1 , CHI?s CHORDNISE POSITIONS OF The FLAP CORNERS, PERCENT of CHORp 

read (R5.5002) Y F 1 » Y F 2 

. ■ . . . YF 1 AND YP 2 ARE DATA AFTER THE RING MAS BEEN YAWED, THE USER 
...♦.MUST DETERMINE THEM FOR EACH SIDESLIP ANGLE, 

READ (R5,5002) CHI1,CHI2 

IF (YFl.GE.YF2) WRITE -(Wh, *201) YF 1 , YF2 

IF (YF1.GF.YF2) STOP 

ETA1R -a f YF 1 - YAVG)/B2 

E T A 1 a £ T A i R 

IF (UNSYH , EQ , 0) ETA* s AMAXi(ETAiR» 0,0) 

FT Ag a ( Yf2 • *AV G)/B2 

ETAIF(l) atTA 1 
ETAIFtaJaETA? 

.....a single flap on a symmetrical wing is allowed to have a 
»•«•* kinked hinge line, this is the reason for having both etai and 
.....etair, if there is a single flap on a symmetrical wing, chii 

♦ ••••APPLIES aT The CENTERLINE, NOT aT E T A l R AS IN ThE OTHER CASES 
PROGRAM ONLY WORKS with THE P A RT OF THE FLAP ON ETa f,T n 
NFL A p S * 1 ... 

IF (UNSYm . EQ. 0 .AND, ETA1R ,G£ , 0,0) NFi ApSs? 

112 CONTINUE 

JJHAX a MAXIMUM NUMBER OF INTEGRATION STaTIOnS 
PP a NUMBER OF CHORDWISE CONTROL POINTS 

» NUMBER OF SPANWtSE CONTROL POINTS ON ENTIRE SPAN 
MN a A REFERENCE NUMBER FQR COMPUTING SPANWISE POSITIONS 

LL a NUMBER OF TIMES THE SPaNWISE INTEGRATION WILL BE DONE 
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c 

READ (R5,50Q3) JJMaX 
READ (R5,5003) PP,nF,NN,LL 

IF (PP.LE.O.OR.NF.LE.O} WRITE . (»6, 6202). PP, NP 
IF (PP.LE.O.OH.NF.LE.O) STOP 
IFfSWTYPE .EQ.OI MM s NF 
IF (MM ,LE. 0} WR I TE ( *6 * 6206 3 
IF (NN.lE.O)STOP 

. .... JJMaX MUST BE am ODD NUMBER 

IF (M0D(JJMAX,2) ,EQ . 0) JJMAX a 2*JJMAX + t 
JJMAX a MIN0(JJMAX, JJMXMX) 
jF ( JJMAX » LE » 0) JJMAX a JJMXMX 

,,,,,'aT THIS POINT JJMaX is the maximum NUMBER OF SPaNWISE integration 
* • * » * T 0 BE ALLOWED 

, JJMAX AND NN MUST 0E RELATED by 
JJMAX+1 b (NN + 1)*2**k 
.....WHERE K s 0 , 1 , 2 ,.,. 

. .... THE FOLLOWING CODE MAKES THE RElaTiONSHjP ABOVE HQLD FOR THE 

.....largest k such that jjm a x is no larger than its current value 

10UH a 2*(NN+1) 

JJMAX1 * JJMAX + j 
JJMAX B NN 
DO 350 J s l, 20 
LLRb J 

IF (IDUM , GT , JJMAX1I GO TO 380 
JJMAX s IOUM - \ 

IDUM s I D U M * 2 

350 CONTINUE 

c ..*,. control should never get to here 

WRITf (Wfe,62<?5) 

STOP 

380 CONTINUE 

LOG 1 a MOO (JJMAX #21 , EQ. 0 

IF (L0G1) WRITE (*6,6203 ) JJMAX, NN 

IF (LOG!) STOP 



iFfLL , IE , 0) LL = UR 
LL » «IN0fLU,iO 

C 

c..., .COMPUTING XSI(LE), CHORDS, CHIFLAP, and C? at A 
C ®E t of INTEGRATION POINTS DfFINfD by jjmax, 

c 

ETAFMX b AMAXl (ABS(ETAIR) , ABS (ETA2) ) 

JJMI * fJJMAX*i)/2 

JJH2«(jjMAx + n/2 

JJMAX1 » JJMAX . I 

NRATlO*CJJHAX+j)/(NN+i) 

C 

DTH£TA*PIO/DFLOAT (JJMAX+1) 

DO 900 J«J, JJH2 
j2xJJMAX-*J*t 

THETAD 6 DTHETA*DFL OAT f J ) 

FT A ( J) B OCOS{THETAD) 

STHETACJJeoSINCTHfcTAo) 

FTAJ«ETA(J) 

ETA(J2)b.ETA(J) 

STheTA ( J2) = STHpTA f J) 

C 

IF (NFLAPS,NE,0) C2 I P ( J ) sC2 I (ETAJ, ETAFMX) 

IF (NFUPS.NE.o) C2IPCJ23 *C2lP(J) 

900 CONTINUE 
C 
C 

C CHICP a CHOROWISE CONTROL POINT POSITIONS (PP VALUES) 

C 

IF(CWTYPE) 1 0R0, 1035, loas 

c 

1035 CONTINUE 
C 

C,,. ..COMPUTING THE CHICP»S ACCORDING TO A SPECIFIED RULE FOR DR WAGNER 
C.....N-INPUT PROGRAM 
C 

PP1 s PP * 1 

IF (? - PP1) 103b, 1036, 1038 

1036 OTHETA a PlD/DFLOATrPP^ 
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00 1037 P ■ 2, PP1 

XSI S «DC0S ( OpLOA T ( P * n* DTHfTA) 

1037 CHICP(P) a Cl, + XSt)/'2, 

J03B C^ICP(I) « 0. 

CHICP(PP) * 1, 

GO TO 1050 
1Q«0 CONTINUE 

,,,, .ARBITRARY CHOHDWISE POSITIONS 

READ(RS#500i) cCHlCPfP)^ P * 1» PP) 

GO TO 1050 
I0«5 CONTINUE 

, . , , ,HULTHOPP f HSU CHOPD^ISE POSITIONS 

IDUM * 2*PP + 1 

DTHETA s PI0/DFL0ATCIDUM) 

DO 1G<*7 I ■ 2 , IOUM, 2 

p a \n 

XSI 8 #0C0S(0THETA*0PL0ATU)) 
io<i7 chicP(P) 8 (I. ♦ xsn/2. 

1050 CONTINUE 

NINDEX 8 INDEX POP SPaNWISE CONTROL POTNT POSITIONS 
IF (SWTYPE ,NE t 0) GO TO 1070 
1060 continue 


...'•.the HULTHOPP SPaNWISE POSITIONS 

DO 1065 N * 1. NP 
1065 NlNOEX(N) b N 
GO TO 1060 
1070 CONTINUE 

.... .USER-SUPPLIED INDICES FOR SPiNWjSE POSITIONS 



IF (IJNSYM.NE.O) NFP b NF 
IF (Ut«SYM,eQ f 0) NFP s ( N F ♦ 1 ) / ? 

READ(R5f5003) (N I NdEX ( N ) , N s 1, NF P } 

IF (UMSYM.ME.O) GO TO 1«P0 
C IF (UNSYM.EU.01 CONTINUE 

DC 1071 I s 1, NFP 

N I M S N F ♦ 1 » I 

NINOEX(NIN) S NF ' ♦ 1 • N INDEX (I ) 

1071 CONTINUE 
1090 CONTINUE 
C 

C CALCULATION OF LEADING EDGE POINTS 4 T INTEGRATION STaTIONS 

C 

IFfLFTYPE.EO.Ol GO TO 117 

CALL COOIM(XLaETaLI#NL * ET A , XLE , JJMaX, XCLE ) 

CALL COOlM(XL»FTALIf NL ,fTaZ,XLZ,1,XCLF) 

GO TO 118 
117 CONTINUE 

CALL ST«aIT(XL,ETaLI»*L , ET a , XLE , JJHAX , ) 

C 

CALL STRaIT(XL»ETaLI^NL » ETaZ , X^Z , 1 , wbi 
IIP CONTINUE 

^ C 

•- C CALCULATION OF TRAILING EDGE POINTS at INTEGRATION STATIONS 

c 

IFCTETYPF.EG.O) go to up 

Call CO0IN(XT, ETATlrNT t ET A t XTE t JJM aX , XCTE 3 
CALL COD IH (XT, fTaTI, NT , FTAZ, XtZ, 1,XCTf 5 
x GO TO I2n 

IIP CONTINUE 

CALL STRaIT(XT,ETaT 1 ,nT ,ETa, XTE, JJMAX,**) 

C 

call strait(Xt,etati ,nt ,etaz,xtz, i ,* 6 i 

120 CONTINUE 

DO 121 msI, JJmAX 
XSILIP( M )sXLE(N)/8? 

C0RDlP(N)8(XTEfM)-XLE(Mn/82 
CORDSG(M) bCOPDjP (M)**2 

121 CONTINUE 



c 

AREAWP b area* 

C8ARWP s CBAR* 

IF CA«eAW.NE, 0,0, AND., CBaRW.NE.Q, 0) GO TO 123 

c 

c DISCONTINUITIES FOR SUBROUTINE INTGWT 

C 

NO b NDL ♦ NOT 

00 1211 I s 1, NOL 
XDlSfd « (YdLU)"YAVG)/B2 
121! CONTINUE , 

DO 1212 I = 1 t NOT 
N s NDt ♦ I 

XDISfN) b (YDT( I )"YAVG)/fl2 
1212 CONTINUE 

IF (AREan,nE,0) GO TO 122 
C 

C COMPUTE AHEAW IF IT IS NOT SUPPLIED 

c , 

WRITE ( W6* 600 l ) 

WRITE (W6»6Q02) 

WR I Tf (W6*6162) 

< « CALL INTGRT tJJMI, jjMAX,ETA r STHETA,CrtRDlP,ND,XOIS f VDIS^6, OU- 

1T2»NLEfT, VALUE) 

A R E A W P a B 2 * * 2 * V A L U E 
C 

CALL INTGRT(JJM 6 X, JJMAX, JJM A X, ETa, STHETa, cordip, NL£FT,YDIS,- 
1 20 1 S , Wh,oUT» T * VALUE V 
AREAWsVALUE*B2**2 
122 CONTINUE 

IF (CBARW,NE t 0,0) GO TO 123 
C 

c COMPUTE CBaRw if IT is NOT SUPPLIED 

C 

WRITE (W6» 6002) 

WRITE (W6#fei60) 

CALL INTr,RT( JJMI, JJMAX, JJKAX^TA,STHETA f COROSO,NO,XOlS, YDIS, W6 f 
10UT2 ,NLEfT, VALUE) 

CBARwP a B**3/(e,0*AREAW )*VaLUE 



CALL INT6RT( JJMAX, JJHaX, JJMaX, ETa, STHETa, CORDSG, NLEFT,Yol$,- 
1 ZOlSf *6, OUT, if VALUE ) 

CBAPW e B**3/(8,0*A«EAW) *VALUE 
C 

123 CONTINUE 

fbref***; *bref**2 

AP s FBRfF 2/AWE AW 
ARP» FRPEF2/ AREARP 
DlFAWeABSCAREAW-AHEAWP5 
DIFARsABS(AR-ARp) 

DIFCBsABS(CBARW-CBAHIwP) 

CBARBP b CBAPw/RrEF 
COROT I eXTfl)-XL(l) 

cordss«xtz w xlz 

C 

THsCOROTI /COROSS 

SSCBARaBREF/CBARW 

DIFSS«ABS(SSCBA»-BREF/CBARwP1 

CBARCDsCBARW/COROSS 

DIFGC#ABSfCBARCD*CBARi«P/CORDSS) 

IF (NFLaPS.EQ.o) GO TO US 

e 

2 C COMPUTE FLAP INFORMATION 

r 

IF(LFTYPE,eG,0) GO TO 12« 

CALL COOIM(xl,ETaLI* nl» ETaIF, XIF, ?,xrLE) 
go TO 125 

124 CONTINUE 

CALL STRaJT (XL# ETaLI, NL, ETaIF, XLF, 2,*M 
C 

125 CONTINUE 
I'P(TETYPE,f0 t O) GO TO 126 

call codiMt, et a ti, nt, etaif, xtf, 2,x c tej 

GO TO 127 

126 CONTINUE 

call straitixt, etati, nt, etaif, xtf, p,w6> 


127 CONTINUE 



non 


CYF(1)'«XTP(1)»XLF{1) 

CYF( ? ) a XTF(2).XlP(?) 

XHH*XLFO)+CVF( 1 )*CHIJ 
XHL2sXUF(2)*CYF(2)*CHI? 

XST1 * XHL1/B? 

XS 12 • XH L 2/B? 

DO 132 M*|,JJMAX 

IF (E7AfM)-ETAlF(2J-i«0E-5) 126*126*131 

128 CONTINUE 

IF<ETA<M)*ETAiF(l)+l # 0E*5) 129,130,130 

129 CONTINUE , 

CORDF(M)»0,0 
CFSG(M) = 0,0 

CHIFPI(M)bCHI1 

GO TO 132 

130 CONTINUE ... 

ET AM e FTA(m) 

XHL«LINIT(fTAH ,XHtl,XHL2*ETAl f fTA2) 

C0RDF(M),XTE(M)-XHL 

CFSO(M)sCORDF(M)**2 

C^lFPI(M)B(XHL-XU(Mn/cO R 0lP(M)/B2 

use average values for flap corners within an epsilon of eta(H} value 

IF (DABScETA(M) . ET A 1 ) , LT , 1 « 0E»5) GO TO 1303 
IF (DABS (ETa ( m ) „ ETA21,LT f i f 0E,5) GO TO 1303 
GO TO 1 32 
1303 CONTINUE 

COHDF(M)sCORDF(M)/2,0 
CFSQ(M)sCFSG(M)/2,0 
GO TO 1 32 

131 CONTINUE 
CORDF(M)sO.O 
CFSQ(M) b 0,0 
CHIFPI (M)b C HI2 

132 CONTINUE 

IF (UN$YM,NE,01 GO TO 1325 
c IF (UNSYM.EU.Ol CONTINUE 

00 1 32U I * 1, JJM2 



n o r> n o o r* 


M = JJMAX + 1-1 

CHlFPIfMJe CHIFPIrl) 

1 32a CONTINUE 
1325 CONTINUE 
C 

IF (AREAF.NE.O.O.ANO.CBARF.NE.O.O) GO TO 13 a 

DISCONTINUITIES FOR SUBROUTINE INTGRT 

DO 1321 1=1, NOT 

XDISC I) = CYOTf I).YAVG)/B? 

1321 CONTINUE 

IF ( NF LAPS, E 0 , 2 ) GO To 1322 
NO = not ♦ ? 

XDIS(NOT + D •= FT A 1 
XOlSfNoT + 2). = F T A 2 

GO TO 1323 

1322 CONTINUE 

NO * NOT + tx 

XDIS(NOTTI) a ETA1 
X0IS(N0T*2) b ETA2 
XOISfNDT + 3) = ■£ T A | 

XOiS (NOT+a) s -EtA? 

1323 CONTINUE 

COMPUTE AREA OF THfc FLaP IF NOT GIVEN 

FOR A SYMMETRIC MING ONLY TH£ A R£a ON ETa,GT ,0 IS CONSIDERED# 

So a Factor of 2 1 5 needed if tuEre is only a single flap 

ARE AFP e AREAF 

IF f AREAF f NE.O) GO TO 133 

WRITE fWS#6002) 

WRIT? (Wfe#M 6 fe) 

CALI INTGRT CJJMI ,JvJMAX,JJMAX,ETA,STHEU f C0RDF,ND,XDlS,YOlS,Wfe. 

1 0UT2, MEfT, VALUE) ' ’ ' * 

AREAFP s 92**2*VA|.UE 

IF (UNSYM.EG.o, AND, NFLAPS,EO, \ \ AREAFP s 2 » * A W£ aF p 

Call INTGRT (JJM&X, J JH A X # JJHAX, ETA# SThETa, COKDF, NLEFT, Ynis, 



1ZDIS, W6,0UT, IND> VALUE) 

AREAF*VA l UE*B2**2 

IF (UNSYM f EQ t O, ANO.NFL aPS.EQ, 1) AREAP s 2,0*AREAF 
C 

133 CONTINUE 

DIFAF s ABS ( ARE AF - AWEaFP) 

C 

C COMPUTE FLAP CBAR IF IT 15 NOT SUPPLIED 
C 

CBARFP a CBARF 

IF CC8ARF ,NE p O) GO TO J3« 

WRITF (Wfc#6002) 

WRITE (W6r6l66) 

CALL INTGFT(JJMI , JJM*X, JJMaX,ETA,STHETa,CFS0,ND,X0IS*Y0IS,W6, 

1 OUT? » NLfcFT, VALUF) 

CBARFP e 8**3/(8,0*AREAF )*VaLUE 

IF {UNSYM, EG . 0 1 ANO.NFL APS ,EQ , 1 ) CBARFP s 2*0* CBARFP 
C 

CALL INTGRT fJJMAK, JJMaX, JJMAX# ETa, STHETa, CFSQ, NIEFT, YOIS, 
1 Z D I S » W6,GUT, T n, D> VALUE) 
cbarf s 8**3/(8,o*areaf )* value 

IF < , UNSYM,EG.o » AND.NFLaPS.EO, , ) CBARF = 2.0*CBARF 

c 

^ 13a CONTINUE 

DIFCBF * ABS (CBARF - CBARFP) 

AfAWsAREAF/APEAW 

C8FCB*bCBARF/CBARW 

DIFCBF s ASSCCBARF - CBARFP) 

DIFAf b ABSfAREAp- * AREAfP) 

LAMDAC b ATAN? f ( X S 1 2 - XSIl)# f E T A2 - ETAl)) 

LAMDCDbLAMDAC*R2D 

C 

C STORE FLAP CORNERS IN ETaF(2»6) 

C 

etaf M>-n«rT*iR 

ETAFfg, i)sETA2 
FTAFfl ,2 )b-ETA2 
ETAF(2,2)8-ETA1R 
1 35 CONTINUE 
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COMPUTE TaNGENTS OF SWEEP ANGLES OF LEADING AND TRAILING EDGES 

DYLMAX s XSlLIPfJJMAX) - XLJNU/R2 

DXUMAX b ETA(JJMAX) ♦ 1,000 

OYLl « XSILIP(I) - XlCn/R2 

OXLi 8 1,000 - ETA ( 1 ) 

DYTMAX 8 "XT (NTJ/B? ♦ XSILIPf JJMAX) ♦ COROIP f JJM aX) 

0YT1 s XSlilPflJ ♦ CORDiP(l) - XT f 1 ) /B2 

LAMtEL(JJMAX) s DYLMAX/DXLMAX 
UAMLERfl) s 0YU/DXU1 
, LAMTEL(JJMAX) a DYTMAX/DX^MAX- 
LAMTER(l) a DYT1/DXL1 

DO 136 Mb1,JJM4XI 

OYLlL * XSILlPfP) • XSjLlPfKM) 

OXILL a ETA(M) • ETA(M+l) 

PTLTL b XSlLlPfM) ♦ COPDlPfM) - XSIilPCMAl) - C0R0 IP(m+1> 

UAMLEI(M)s OYULL/DXULL 
LAMTELfM)* DVUTt/DXLLL 
L AM^FK C M ♦ 1 ) b "L Am LEL( m 5 
LAMTERfM*!)* •LANTFL (M) 

136 CONTINUE 

NC e NUMBER OF CHORpwiSE POINTS AT WHICH PRESSURES WILL RE COMPUTED 

IFfCPRESS ,EQ t 0.OP,NFLAPS,EQ,0) HO TO 1700 
READ (R5,5003) NC 
IF ( NC ,IT • 1 «DR,Nr »GT ,50 ) N" s <19 
C 

C, .COMPUTING THE PRESSURE COEFFICIENT for THE FLAP 
C,,,., M ODF OF PRESSURE AT THE SPaNMSE CONTROL POINT LOCATIONS 

C 

WRITE {WferSOOl) 

WRITF (W6r600«) 

OTHpTA b PI/FLOAT (NC+1 ) 

00 1520 N 8 1, NC 



THETA x FLOAT£N)*DTHETA 
X(N) s *COS(THETA) 

CHI C N ) x (y tN) + l,)/2, 

1520 S Q x f m J x S I Ki f THfcT A) 

WRJTE . fW6#6035) ( CM 1 ( N) , Nx l , NC ) 

WRITf £W6»6Q36) (X (N) , tUl , NC ) 

PAIR x NFL APS ,GT, 1 
C08UC x COS(LAMDaC) 

C 0 N 1 x SQRTM ,-CM 4 CH*cOSlO**2)/COSLC 
CON2x-2,/(PI*COni ) 

SlGNxi, , 

WRITE (W6,6Q01) 

WRITE (W6# 6004) 

WRITE (W6f6036) 

1550 CONTINUE 

DO 1600 J * NRATIO, JJHAX , NRaTIO 
FTAJ x ETA(J) 

WRITE (W6» 6037) ETAJ 
OURV x 2,/COR0IP(J) *rOMi 
YlR x DUMY*(ETA.t - ET A 1 R ) 

Y2«. « DUMV*(ETAJ - ETA2) 

Y1L x *DUMY*(ETAJ+ETA1R) 

Y2t « •0UMV*(ETAj*ETA2) 

XC x 2,*CHIFRI (J)- l, 

3QRTXC x SoRT ( \ *«XC*XC) 

DUMY * C0N2*C2tP(J)/SGRTXC 
00 1560 Nx 1; NC 
XNXC * X(N)- XC 

IF(ABS<XNXC).LT.i .E-7) XNXC*i,E- 7 
TIP x 0(XNXC# -Y1R) 

T2R x Q(XNXC* -Y2R) 

Tit « QfXNXCt ^ V 1 L ) 

T2i x QCXNXC, -Y2|_) 

DCP(N) x D U N Y * SQX(N) * ALOG( T1R/T2R ) 

1560 IF(PAIR) DCP(N) x DC P C N ) > DUmYaSGX < N) * AlOG ( T H/T 2 L) *5 If,N 
WRITE (W6,6041) (DCP(N)#N«l f NC) . 

1600 CONTINUE 

IF (.NOT.PaIR .OR, SIGN.LT, 0.) GO TO 1700 
WRITE (W6t600i) 



WRITE (Wfe,600<O 
WRITF (W6, 6039) 

SIGN a »i f 
GO TO 1SSO 
1700 CONTINUE 

C 

C LINE PRINTER OUTPUT 

C 

WRITE (W6,6001) 

WRITE (W6,6004) 

WRITE (W6»6005) , 

WRITE ( W6# 6 020 ) 

WRITE (W6> 6009) AREAW, DIFaw 
WRITE(W 6,630i) BRER 
WRIT E (W6f 6006) B2 
BRAT«BREE/B2 ' 

WRITE (W6> 6302) BRaT 
WRITF (W6f 60 1 1 ) CBaRW, OlfrCB 
WRITE (W6 f 6010) CORDSS 
WRITE ( W6# 60 1 2) A R # DIFAR 
WRITE ( w6# 60 1 3 ) TR 
WRITE (W6f601tt) SSCBAR, DIFSS 
WRITE (*6,6015) CBaRCD, DIFGC 
WRITE (W6,6l5l ) MaCH 
WRITE (W6#6002) 

WRITE (W6»60U2) 

WRITE (*6,6001) 

WRITE (W6,6Q02) 

WRITF (*6,6016) LFSCAL 

WRITE (*6,6017) (Yim,Xl(I),I»l,NL ) 

WRITE (W6,600g) 

WRITE (*6,6018) TESCal 
WRITF (W6,6017) ( YT ( I) , XT ( I ) , I a 1 , NT ) 
IF (NDL.EO.O) GO TQ U01 
IF (LETYPE.f-O.O) GO Tn 1601 
WRITE (W6,6001) 

WRITE (*6,6002) 

WRITf (*6,6170) 

WRITE (W6,6G02) 



nr*c% r*no r>r*n 


WRITE (W6,M72) (YOL(I)r IbI'NDL) 

1801 CONTINUE 

IE (NOT.EO.O) GO TO 1802 
IF (TETYPE.fcQtO) CO TO 1802 
IF (LETYPE.EO.OJWRITE (Nt'tOOl) 

WRITE (Wb,6002) 

WR1TF (Hbtbni) 

WRITE (W«»,b002) 

WRITE (W6 f fei72) (YDT(I) # I«1,NDT) 

1802 CONTINUE 

if (trnang.eo.o.i go ;ro i ac>3 

WRITE (Wfe,bO<li) 

IF (TRNfRM.LT.o) WRlT£CN6»hftR<0 T«N*NG 
jF(TRNFRH,FQ.O) WRfTE ( Wfe, 60R5) TRnANg 
IF ( TRNfRM, 6T,01 WRITE (Wb,60afc) TRNANG 
18(53 CONTINUE 

IF CNFUPS'EQyO) 50 To i ?a9 

PRINT THE FLaP DATA 

WRITE -rW6#6001) 

WRITE CW6#fc00<O 
WRITp (*6,6112) 

WRITE (Wb,6i08) tAHDCD 

WRITE (N 6,-6 114) aFaw, C8FC0W 

WRITE (W6,6135) AREAF, DIFaF, CBARE, DTFCbF 

IF(Nfl,AP8,Ftt. 2) GO TO 1 247 

If (UNSYM.fO.O) GO TO 12ab 

UNSYMMETRICaL wing - SINGLE FLAP 

WRITE (W6,616n 

WRITE (WfeifeHS) E T A i R » ETA2» CHIl, CHI2 
WRITE (*6,6117) CYrU), rYF(2), XSI1, XSI2 
GO TO 12U9 
1246 CONTINUE 


single symetrical flap 
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WRITE (*6,6132) ETAlR, ETA1, F TA 2 , CHI2 f CHI1, CHI2 
WRITE (*6,6133) XSI2, XSI1» XSI2, CYF(?), CYF(l), CYF(2J 
GO TO 1249 
J 24 7 CONTINUE 

SYMMETRICAL WING • TWO FLAPS 
WRITE ( W6# 6 1 29 ) 

WRITE (*6,6115) E T A 1 » * ETA2, CHJi, CHI? 

WRITE (*6,61 17) CYf(1)» CYF (2) » XSII, XS 1 2 
WRITE (*6, 6131 ) , 

WRITE (*6,6115) ETaF(i # 2), ETaF( 2,2), CMl2, CHlj 
WRITE (*6,6117) CYe(2>» C Ye ( 1 ) , XS 12 , XSII 

1249 CONTINUE ' 

IFtlPlNEO ,EO. 0) GO TO 1255 
WRITE (W6, 6.00 1 ) 

WRITE (*6,6004) 

IE (UNSYM,fQ t Oj JJMP x JJM2 
IF (UNSYM.NE.O) JJMP „ JJMAX 
IF (NFLAPS.EQ.O) GO TO 1251 

WRITE (*#,,6107) ' 

DO I25fl J«l,jjHP 
ETAJsETA(J) 

WRITE (W6,6105) ETaJ, XSILIP(J), COROIP(J), C.HIFPlfJ), C?I p (J), 

1 LAMLfcL(.j)» LAMLER ( J) , l A m TEL ( 1) , LAMTER(J) 

1250 CONTINUE 
GO TO 1255 

1251 CONTINUE 
WRITE ( *6 , 6 1 58 ) 

DO 1252 J s 1, JJMP 
ETAJxETA(J) 

WRITE (W6,6105) ETaJ, XSILIP(J), CORDIP(J), LaMLEL(J)# LaML 

1 E R f J ) , LAMTEL (j) , LAMTEKfj) 

1252 CONTINUE 
1255 CONTINUE 

WRITE (*6,6001) 

WRlTg (*6,6002) 

WRITE (Wfe,6109) 

00 1260 Psl,PP 
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WRITE (W6,M05) CHICP(PJ 
I2fe0 CONTINUE 

WRITE ( W6» 6004) 

WRITE ?Wfc,feliO) 

WRITE (Wfe»6002) 

DO 1300 N*l, NF 
INDEX*NINDEXfN)*NRATlO 
ETAN 9 ETACINOEX) 

WRITE (WR»fel05) CT AN 
1300 CONTINUE 

WRITE(W6,6305) JJMaX , 

UNIT 7 OUTPUT 

IF (N07 ,NF, 0) WRITf (*6,6303) 

IF CN07 ,NF, o) GO TO 2j3 

» * • » « GEMEIt, IS A SUBROUTINE WHICH COMPUTES THE GEOMETRY FILE NaME AND 

ISSUES aPPPORRIaTE DDEF COMMANDS, THIS SUBROUTINE IS TO hE USED ONLY 
, , • t » 0 N THE AMES’ TSS COMMUTER SYSTEM, 

CALL GEMFIL(ID) 

REWIND NT 

IF (NDL ,LE, o) GO TO aia 
00 211 I * 1, NDL - 
YDLfl) « YDL(I)/b? 

all CONTINUE 
2iR N0(_ « MAxO (NDl* 1 ) 

• IF (NOT ,LE, 0) GO TO 215 
DO 212 I * 1, NOT 
YDT ( I ) t YDT f t ) /B2 

2 l 2 CONTINUE 
2l5 NOT * MAXQfNDT, 1) 

WftlTg ( W 7 ) ID* PP, Nf# CWTYPf, SWTYPf, UNSYM, NOL * NOT, NN, JJMAX, - 
1 NfLAPS, TITLE, NTITL 

CbCBARBR 
BRbBRAT 
k*nf 

p*pp 



nnors 


NsNDL 
l«NDT 
J s J J M A X 

^RITE(M7) (CHICP(X) > X«1 f P) r (NINDEX(I) »!■! f K)f f LAMIEL ( I ) , 1=1 » J) , (tAM 
tLERtDf I«UJ)» (LAMTEL(t) # XBl r J), CLAMTER ( I) , 1 = 1 , J) , 

i(FTAa)»isi#j), ( stmeta(i) , tsi , jj , <xsilzpcz)#t s i#J)# ccordi 

lP(I)»I = lfJ)fBR#C»AR,TH#MACH,(yDL(l),!cl,Ki),(YDT(I),rsl,t) 
ip (mplaps,me,o)write(N7) mach, lamdac, cbfcbw, afaw, e t a i * etas, 

1XUF, CYf,XSI 1, XSJ2, (CHlpPIfI)#I*l,JjMAX), (C2IP ( I ) , jsl 

lfJJMAX) 

ENOFILE W7 
WRI TE (H6, 6 3 0 <4 ) 

213 CONTINUE 

dounwash values fur unit 8 

alfab b array of oownhash values 

WRITE (Wfr,600n 

NS V M 3 0 

NASYM 8 (5 

DO 179 I B 1, 10 

IP (BCS(in NS YM S NS VM * 1 
if cecAScm nasym « nasym + 1 
179 CONTINUE 

IF (NSYM.EQ.O, AND.NASYM.EO.O) WRITE (W6,6174) 

IF (NSYM t EQ,O t ANO,NASY*,EQ.O) GO TO 19ft 

LOG1 b N06 ,EQ, 0 
ID3b0 

MMp*(NFM)/2 

iFCUNSYHjNE.OjHMPaMM 
MMP a*NF/2 

IF(UNSyH,NE,0)MMPA*Mm 
DO ISO 

INOeXbNINDfX(M) 

ETACP(M)«Ef A(INOEX) 

180 CONTINUE 

iFf.NOT.LOGl) GO TO 1809 
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. . BCPlt IS A SUBROUTINE WHICH COMPUTES THE BOUNDARY CONDITION FI^E NAME 

and issues appropriate ddef commands; this subroutine is TO be USED only 

*||"on THE AMES' TSS COMPUTER SYSTEM. 

c 

103*0 

. CALL BCFILC 101,1035 

IF (LOG'S) REWIND WB . . . 

IP ( L. 0 G 1 ) WRITE (W8) ID»I03»TITLE. UNSYM, NSYM, NaSYM, BCS, BCaS, 

1 PP* CWTYPE ^ NP> MMP. MMPa f SWT YPE f (C hICP(P) t P«t t PP ) t 
2 ( NJ NOEX ( M) »M»l#MMP)f(ETACP(M) f M»l f MMP) 

C i 

1809 CONTINUE 

WRITE (W6# 6002) 

IF (LOGI)RRITE (*6,6175) 

IF ( ,N0T,L051)wrITE (W6,6300) 

WRITF ( w<> f 6.002 ) 

.. nmax C NF 

' " v IF (UNSYM, EG. 0) nmax b ( NF ♦ i ) /2 
11 C .NOT ,BCS ( 1 ) ) GO TO 182 
C 

C UNIFORM DOWNWaSH MODE 

® c 

DO 181 N * 1. NMAX 
00 181 P « 1. PP 

A L F A B ( P f N ) 8 1,0 

181 CONTINUE 

IF (LOGllWRITE ( W8 ) ( ( ALFAB(P»M »P=1 , PP) , Ns * , NmAX ) 

WRITE (W6,6lO«) 

WfilTF (W6,6152) 

WRITF rW6,61Q4) 

ASSIGN i82 TO LOBO 
CO TO 301 

182 CONTINUE 

IF ( t NOT t BCS(2) ) GO TO 164 
C 

C PITCHING MODE 

C 

DO 183 N s 1, NMAX 
00 183 P * 1, p P 



JKA s NRATI0*NJNDEX(N} 

ALFAB(P# N) • (XSILIP(JRA) ♦ CMICP(P)*rORDIPf JR a) )*SS C 8ar 
183 CONTINUE 

If (L06l ) WRITE (we) ((ALFAB(P*N)fPa|»PPWNal«NMAX) 

WRITE 1*6,6104) 

WRITE (W6#fcl53) 

WRITE (W6,6104) 

ASSIGN ;8a TO LOftO 
GO TO 301 

lea continue 

If t.NOT.BCSd)) GO TO 186 
C 

C LINEAR SYMMETRIC TWIST 

c 

DO 185 N s l, NMAX 

00 185 P # i, PP 

JR A s NRATlO*NINDFX (N) 

AIFAB fP,N) *DA8S (ET A ( JRA ) ) 

185 CONTINUE 

If (LOGl ) WRITE ( W'8 ) ( (ALFAB(P,N) ,P*i ,PP) ,Nei ,N mAX) 

WfTITF {*6,6104) 

WRlTF fW6,6157) 

WRITE f W6,6104) 

ASSIGN | 86 TO lOpO 
GO TO 301 

186 CONTINUE 

IF {,NOT t BCS(«)) GO TO 188 
C 

C CALL TO FLPOWN 

C 

CALL f LPDWN(CWTYPE, LL, Pp, NF, Mach, NFLaPS, NN, JJMaX,TITL» 

IF, LAMDAC, NjnDEX, CwICP, LAMlEu LAMlER, EAMTFl, LAMTER,ETaf« 

1R5, W6,ETA, 

1 STmETa, HSILIP# CORDIP# CHIfPI , C3IP, ALPHA) 

IF (NFLAPS.EQ.l) GO TO 1872 
DO 1 6 7 i m * t, NmA x 
NfN b Nf 4 1 • N 

DO 1871 P B 1, Pp 

ALFAB(P,N) a m( ALPHA (LL ,P , N ) * ALPH A (LL , P , NF N ) ) 



non 


1071 CONTINUE 
GO TO 1 87 

1072 CONTINUE 

IF (UNSYM.fQ.O! GO TO 107a 
00 1873 N » If NMAX 
DP 1873 Pel, PP 
ALFAB(P,N) «*ALPHA(U. f P,N) 

1073 CONTINUE 
GO TO 107 

1874 CONTINUE 

DO 1875 Ns j, NmAX 
DO 1075 P«1 , PP 
ALFAB(P,N) s-ALPNA (U # P, N) 

1075 CONTINUE 
187 CONTINUE 

Nft I TF (Wfefblfto) 

WRITf t*6ffcl5«) 

WRITE (W6,fcl0«) 

ASSIGN 10 B TO uoco 
GO TO 2ot 
■ 188 CONTINUE 

aLFab(Pp,nm a x) read from cards 

DO 109 I * 5, 10 
IF ( f NOT a BC8'(I)) GO TO 169 

READ (B5/5002) ((ALFABCP»N) f Psi»PP)»Nsi»NMAX) 

IF (LOGUNRlTp C«6) ((ALFAB(P»N)fPel»PP)rNslfNMAXl 
WRITE (W6 t.61.04) 

WRITf (*6»6i56) 

WRITE 

Do 1691 N s 1 , Nm A V 

WRITE fW6#6105) (AUFAB(P,N) ,Pel ,PP) 

1891 CONTINUE 
109 CONTINUE 
190 CONTINUE 

IF (UNSYN.FU.O) nma * * *P/2 
IF (,N0T,8CAS( 1 ) ) GO TO 192 
C 
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C UNIFORM ANTISYMMETRIC MODE 
C 

DO 191 N * 1, NM AX 
DO 191 P s i, PR 
JRA s NRATI0*NJNDEX(N) 

ALFAB(P,N) s DSICN(1 t OO,tTA(JRA) ) 

IF (Dabs (ET A (JRA) ) ,LT, 1,0-5) AlFaB (R, N) bo’. 

191 CONTINUE 

IF (LOGl)NRITE (W8) ( ( AlF AB (P , N) , Pal , PP) , Nn 1 , NMAX ) 

WR I TE (W6,6l04) 

WRITE (W6, 6159) 

WRITE (*6,6104) 

ASSIGN 192 TO LORO 
GO TO 501 

192 CONTINUE 

IF (,N0T*.BCAS(2)) 50 TO 194 

LINEAR ANTISYMMETRIC TWIST 

DO 193 N b 1, NMax 
DO l 93 P a i, PP 
JRA » NRATIObNJNOfX (N) 

ALFAB(P,N) r -ETA(JRA) 

193 CONTINUE 

IF (LOG 1 ) WRITE (W8) (CALFAB(P#N),P«i ,PP),N«i,NmAX) 

WRITF (W6,fel04) 
write (W6,6l60) 

WRITE (W6,6104) 

ASSIGN i9« TO LO0O 
GO TO 301 

194 CONTINUE 

IF ( ,N0T,BCAS(3) ) GO TO 196 
CALL TO FLPDWN 

IF ( ■ NOT , BCS(4))CaLL FLPDWNfCWTYPE, LL, PP, NF, MaCH, NFLaP 

is, nn, jjmax, title# lamoac, njndex, chicp.lamlel, lamler, lamtel# 
1 LAMTER, ETaf,R5,W6,ETa, SThETA, XSlLlP, CORDiP, 

1C2IP, ALPHA) 



r% o r> 


DO 195 N * 1* NMAX 
NN fJf ♦ 1 * N 

DO 195 P * 1, PP 

ALFAB(P,N) b w(ALPHA (LL,P,Nl * ALPHA (LL# P# NFN) ) 

195 CONTINUE 
WRITE (W6,6104) 

WRITE (*6,6155) 

WRITE (Wfe,felf)4) 

ASSIGN 196 TO LOCO 
SO TO 201 

196 CONTINUE 

ALPAB(PP,NM*X) READ FROM CAROS 

DO 197 I *4, 10 
IF ( t NOT t BCA$(l)) GO TO J 97 

READ CR5,5002> ((ALFAB(P,N),Pa l# PP),NBi ( NMAX3 
IF (L0G1)WRIT E (WB) C(ALFAB(P^),P«1,PP),N»1.NMAX) 

WRITE (W6, 61041 
WRITE (Wfc,6l5fc) 

WRlTg (W6,61O0) 

DO 1971 N s 1, NMAx 
WRITE ( W6 » 61 05) ( ALF AB ( P , N) , Ps 1 , PP j 
1971 CONTINUE 

197 CONTINUE 
C 

198 CONTINUE 

IF (LOGO ENnFILE we 
STOP 
c 

201 CONTINUE 

00 1400 P « 1, PP 
DO 1400 N b 1, NmAX 

1 s NRATIO * NIND&XCN) 

TDUM s CHICP (P3 ,GT. CHJFPI(I) 

Z * ETA(t) 

IF (TOUm ;aN 0. (ABSfZ-ETAl) ,LE t l.E-5 ,0R,A6S(Z*ETA2) ,LE , 1,E*5)>* 
1 ALFAB(P,N) a ALF*B(P#N) ♦ 0,5 
TOUR b TDUM .AND, SNGL(ETAfl)) ,LT, ETA2 - l.E-5 



n n c* r* r* rs n 


TDUM s TDUM .AND, 8N6U(ETA(D) ,GT, ETA 1 ♦ t,E*5 
IFfTDUM) AlFABfP,*) a ALFAB(P,N) + t',0 
14 0.0 CONTINUE 

IF (LCGl)WRITE (WS) ((ALFAB(P,N),Pal,PP),Nal,NMAX) 

DO 1510 N ■ Nmax 

WRITE (*6,6105) (AuFABfPiN^Pel ,PP) 

1510 CONTINUE 

GO TO LOCO,' 1186,196) 

301 CONTINUE 

00 302 N s 1, NMAX 

WRITE fW6#6l05) ( ALF AB (P, N) , Pa ) , PP ) 

302 CONTINUE 

GO TO LOBO, {182, 104, 186, 190, 192, 194) 

FORMAT statements 


•*,j.TME FOLLOWING FORMAT DEPENDS on the integer word length 
* • • « • OF THE COMPUTER BEING USED', 

« ■ * • , Change t h e following format as required by integer word length 

5001 FORMAT { 2 0 A 4 ) 

5002 FORMAT f 8F 1 0 , 0 ) 

5003 FORMAT ( 1615) 

5000 FORMAT(lOLl) 

5005 FORMAT C2l5, 3FiO,0) 

C 

6001 FORMaT{1H1) 

6002 FORMAT (IMq) 

6003 F0Rm 4 T(1H ) 

6004 Format (////) 

6005 F0RmaT{1X,40HGE0METRICAL data FOR the FLAPPED wing ) 

C,,,,,CHaNGE the following format as rfquireo by INTEGER WORD length 

6006 F0RMAT(26H0TlTLfc OF CONFIGURATION a , 7X, 20A4) 

6007 FORMAT { 33HoaSSIGn£ 0 IDENTIFICATION NUMBER a ,15) 

6008 FORMAT{37H 02* EFFECTIVE SEMI. SPAN ,lHs,F12,5) 

6009 aORMATflX, 12 h wjNg AREA ,24X,1 ms , f 12,5, 3X,lH( » Fl2,5,lH) ) 

6010 F0RmaT( 1X# 13H ROOT CHORD , 23X , 1 Ha , F 1 2 , 5 ) 

6011 FORMAT (37H CBaR, LONG, REF, LENGTH 

l H ( , Fl2.5,lH) ) 


»lHa,Fi2,5,3X, 1 



6012 F 0 R M A T e 1 6 H 0 ASPECT RATIO , 2iX,tHa , Fl?,5, 3X,1H( .F12.5, 1H)) 
60}3 FORMAT ( 1 x » 1 4H Taper Ratio , 22x,lHa , F12.5) 

6014 FORMAT (IX, 33H RREF/CBaR , 3X,lHa , Fl2,b, 

1 3X » 1 H ( , F12.5.1H) ) 

6015 FoRMAT(1X,34H CBAR/RQOT CMq«D , 2X,lHs , Fi2,5, 

1 3x » 1 H f , F12,5, 1H) ) 

6016 Format ( IX# 26H LFADInG edge CO"ORDItv*TES ,SX, 16H(SCALE FACTOR s 
1/ MO. 6# 1 M) // , 1 1 X , 2hYL» UX.2WXI) 

6017 F0RMAT(2F13,5) 

6018 FORMAT ( 1 X f 28H TRAILING EDGE CO.ORD I N A TES , 5 X , } 6 H f SCALE FACTOR a 
1 , Fi0,6, lH)// r 1 IX# 2HYTf 11X,2MXT) 

6019 FoRMATf 1H0/31H0FLAPPED MING GEOMETRY PROGRAM /31H 

1 SSBSSSSSESSSSSS //) 

6020 FORMATf lHO/llHOwING DATA / 11 H --------- /1M0) 

6021 FORMAT (// l oHOC^TYPE a ,13/lQH SRTYPE a ,13/lOH IPINFO a ,13/iOH CP 
1RESS a , 1 3 / 1 0 h UNSYM a ,13/iOH TRNpRM a ,13 /IOh OUT * ,I3/10h 

,, 1 N07 a -,13/lOH N06 a ,13) 

C 

6035 FORMAT ( /«9HOCHORD*I5E STATIONS FOR THE PRESSURE COEFFICIENT 

1 t 31h(PERCEnt LOCAL CHORD FROM L,E, ) ,//, C2P10f12,5) ) 

C 

6036 FORMAT ( /49M0CW0RDWISE STATIONS FOR THE PRESSURE COEFFICIENT 

. -1,R5h(fRACTiQN Of LOCAL SEMI-CHORD pROM miO-ChQRD ) , / ( 1 Of 1 2 ,5 ) ) 

c ‘ 

6037 FORMAT ( 26H0SPANM9F POSITION, ETA a , F8.5) 

C 

6038 FORMAT (IHO.tt 3 H PRESSURE COEFFICIENT FOR THE SYMMETRIC FLAP, 

1 l^H RReSSURf MODE,///) 

c 

6039 FORMAT (1X,49HPRESSURE COEFFICIENT FOR THfc ANT W5YMMETR IC FLAP, 

1 1 R H PRESSURE MODE,///) 

C 

6041 FORMATf 10M2, 5) 

6042 FORMAT { 1 X , 5i HTHE NUMBERS IN PARENTHESES a»E ERROR ESTIMATES 

1 /30w (SEE PROrRAM DOCUMENTATION) ) 

6043 FORMAT (SOHOTHE ABOVE VALUES HERE DETERMINED BY APPLYING ) 

6044 FORMAT ( l X , gH a . F7,5,«9H DEGREE CRANKING TR A NSFORM 4 T ION TO THE IN 

1 PUT ) 

6045 FORMAT ( 1 X , 2HA , F7.3.49H DEGREE SHEAR I mG TRANSFORMATION TO THE In 



iPUT ) 

6096 FORMAT (1X,2HA , F7,3.fl-9H DEGREE ROTATION TRANSFORMATION TO THE IN 
1 PUT ) 

C 

6109 FORMAT (////) 

6105 FORMAT (9F ja,5) 

610? FORMAT (lX,RaMGEOMETRIeAU DATA *T THE INTEGRATION STATIONS, 

1 /jy, 10X, 3hEt A , 7X, 6hX(LE)/R2, 9 X i 5hC/B2, 5x, 15h(x(LE)- 

1X(TL))/C» 5X » 2HC2, 10X, 6hTANLEL, 8X, 6mTANLER,8X, 

1 bHTANTELf 8X, 6MTANTER) 

6108 FORMAT (29HQMINGEL XNE SWEEP * ,F12 t 5, 3X, 7M0EGREE5 ) 

6109 FORMAT (iX,5lHCH0RDWISE CONTROL POINT LOCATIONS REFERENCED TO THE, 

l i 2 m local chord,///) 

6110 FORMAT f lX, 50H SPaNWISE CONTROL POINT LOCATIONS REFERENCED TO B2 ) 

6112 FORMAT C 1 OH FLAP DATA /lOH /JH01 

6115 FORMAT ( 1 X , 37HSPANW ISE LOCATION OF FLAP EDGES/B2 , / , 2F 1 2 , 5 , / , 90H 

1 RELATIVE ChORDWISE LOCATIONS Of CORNERS, /2rl 2,5) 

6117 FORMAT ( 32H CHORD LENGTHS AT FLAP EDGES a«E , / , 2F 1 2 , 5 , / , 

1 25H XSI'S OF THE CORNERS ARE, /, 2f t 2 , 5) 

6129 FORMAT (1X,//,20H FOR FLAP NUMBER ONE) 

6131 FORMAT (iX,//,?0h FOR FLAP NUMBER TWO) 

613? FORMAT (lX,//,a6H DATA FOR THE SjNGlE Flap (L»H f S, , CENTERLINE,, 

1 7h R,H«S) , //,25m SPANWISE P0SITI0NS/B2 , / , If 1 5, 5, // , 

1 52h ChORDWjSE LOCATIONS (PERCENT LOCAL CHORD FROM L.E.), 

1 /,2P3f15,5) 

6133 f 0 rmat (ix,//,37H cHnROwisE locations (actual pqs/&2) /3fis,5, 

1 //,19m CHORD LENGTHS, /,3f15, 5) 

6139 FORMAT (?9H flap AREA/WING AREA e , F12.5/, 24H FLAP CB A 

IR/wING CBAR * , f 1 2 , 5 ) 

6135 FORMAT (IX, 23HFLAP AREA s , F12.5, 3X, l H ( , F 1 2 , 5 , 1 H ) , / 

1, 24m CBAR OF THE FLAP s , f 1 2 , 5 , 3X, 1 H C , Fl2,5, 1 h)) 


6151 FORMAT ( 2 X , 1 2 H MaCH NUMBER , 23X , i H* , F 1 2 , 5 ) 

6152 FQRmaT ( t X, aOHUNjFORM DOWNWaSH MODE ) 

6153 FORMAT (IX, aOHPiTCHjNG MODE ) 

6159 FORMAT (IX, «OHC0MPUTED BY FLPDWN (SYMETRIf MODE) ) 

fe l 55 FORM A T (IX, aOHCOMPUTED PY FLPDWN ( aNT I-SYMETRIC MODE) ) 

6156 FORMAT (IX, 90HDISTRIBUTIONS REaD FROM CARDS ) 

6157 FORMAT (IX, 9QHLINEAR SYMMETRIC TWIST ) 


6 I 58 FORMAT (lX,aflHGEOMETRlCAL DATA aT THE INTEGRATION STATIONS, 



9X, 6HTANLEU. 


1 //, tOX, 3MET A » 7X» 8HX(LEJ/B2, 9X, 4HC/B2, 

18X, 6hTANLER, 8X, 6hT ANTEL t 6X, 6HTANTER) 

6159 FORMAT (IX, aOHONJFORM aNT I -S YMMETR I C ) 

feUe FORMAT (IX, uOHLINEAR ANTI-SYMMETRIC TWIST 5 

6161 FORMAT (IX# //, UOH DATA FOR A SInGl? ASYMMETRICAL FLAP ) 

6162 format (ix,apM integrations for the wing area i 

616a FORMAT dX'ttOH INTEGRATIONS FOR THE WING CRaR i 

6166 format (ix,aoH integrations for the flap area i 

6168 FORMAT (1X,40M INTEGRATIONS FQP THF FLAP CBaR ) 

6170 FORMAT ( jX, uflMTHE WIng HaS KINK$ IN THF. LEADING EDGE , «- 

10HAT The FOLLOWING ETA LOCATIONS ) 

6171 FORMAT (IX,«0HTHE WlNG HaS KINKS IN THE TRAILING EDGE , a. 

iohat the following eta locations i 

6172 FORMAT (FI5.5) 


6174 FORMAT ( 4*HqN0 BOUNDARY CONDITIONS WFRE STORED ON UNIT W* ) 

6175 FORMAT ( 42H0THE FOLLOWING SETS OF BOUNDARY CONDITIONS 

1 40H WERE STORED ON UNjT «8 ) 

C 

6201 FORMAT (IX, 29HBAD FLAP GEOMETRY YPt,GF,YF2, 4X, 6HYF1 = , Fl5,5, - 

l 4X, 6 HYf2 b , F 1 5 , 5 ) 

6202 FORMAT ( 1 x i 42HNUM»ER OF CONTROL POINTS MUST BE PGSlTlyE, 4- 

IX, 6hPP * f 1 5 # 4 X r 6MNF s , J S 5 

g 62o3 FORMAT ( 50H0 ERROR — JJM4y MUST BE MADE AN 000 NUMBER /- 

UOH JJMAX s , 16, 10X, 5hNN e , 16 ) 

6205 FORMAT flX, 8H.ST0P 3501 

6?06 FORMAT (iX^ 541HNN MUST BE ASSIGNED A VALUE WHEN SWTYPfc IS NQT 7f RO - 
1 ) 

6300 FQRmAT { 72H0THE FOLLOWING SETS QF BOUNDARY CONDITIONS WERE COM- 
1PUTFD BUT NOT STORED ) 

6301 FORMAT (37H BREF, LAT. REF, LENGTH ,1H«, Fi2,5) 

63q 2 F0RMATC11H BREF/B2 , 26X, iH» f Fl2 f 5) 

6303 fORMAT(//47hONO GEOMETRY INFORMATION WAS STORfD ON UNIT N7 ) 

6304 F0RMaT(//5«HqTHE GEOMETRY INFORMATION HAS BEEN STORED ON UNIT W? - 
1 ) 

6305 F0RhAT(///IX,7HjjmAX s , 14////////////} 

. end 
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END END END END END END END END END END END END END END END END 


SUBROUTINE CODIH CYI#*I,NI,T, ANS,NA,XK1 

e 

C**** A controlled deviation ITERPOLaTION method 

c 

double PRECISION t 

DIMENSION XI(NI), YI(NI), T(NA>, ANS(NA) 

c 

c 

N*N I 

SIGN * 1,0 

IF (XI(Nt),LT,XXf 1)) SIGN a .1,0 

DO 910 l£* 1 ,NA 

X*T(IE) 

1 00 IFfN-2Jil0,i20,200 
110 V a VJ (N) 
vo GO TO 900 

w 120 v a (Yi(g).vm))/(xi(2).xnn)* (x.xm)) +yim) 

GO TO 900 
200 J = 1 

210 IF ( S 1 GN* ( X I ( J ) - X)) 230,220,2*50 
220 V aYI(J) 

GO TO 900 
230 J ■ J+l 

IF (J"NJ210, 210*250 
250 lF(J-2) 120, 155,260 
155 J a 3 
JJ a } 

GO TO ?85 

260 IF(.J.N)280, 265*270 
265 J * N»1 
JJ a 2 
GO TO ?85 

270 YB (YI(N)-YnN.i))/CXHN)-XI(N.in* CX.Xl(N.imrTfN.l) 



GO TO 900 
2fi0 a i 

285 IF(N.3)29Q,290,295 
290 J s 3 '■ 

295 K a 

m ■ K-t 
L B J* 1 
At s X*XJ(M) 

Ag b X« X J ( K ) 

A3 b X*XI(J) 

AL s (X*XI(K))/(XJ(J)*XHlO) 

s • Au*yicJ)+(i.o*AL)*yi(K) 

ClB A3*A2 /((xi(M)»KI(K))*(XX(M)bXI(J))) 
C2b Ai*A3/((Xl(K)»XX(P))*(XX(K)»XIfJ))) 
C3 b 'A2*Ai/ 4(XX(J)*XHM))«(XX(J)-»XX(K)) ) 
PI s Cl*yi(M)+C2*Y!CK)+C3*Vi(J) 
IF(N*3)305,305,310 
3fl5 P2 b Pi 
GO TO 31 5 
310 A4 b X*X! (U 

C ii* A4*A3/f (XI(K}*XI(J) )*(XI(K)-XI (L) ) ) 
C5b A2*A4/({XifJ)*Xl(K))*<xi<J)wXI(L)>) 
Cfes A3*A?/((XI (L)-XIC*) )*(XI CL) -XI f J) ) ) 
P2 s C4*YI(IO*C5 *Yt(J)+CMYICU 
3 1 S GO TO f320,330,350), JJ 
320 p 2 » p t 

AL s tX*XI(l))/<XI<2)«XI(l)) 

S ■ aL*YI( 2)+ (l.«M)*yi(i) 

Pi* s + XK*(P2-S) 

GO TO 350 
330 p l * P2 

AL * CX-XICNiil))/CXI(NipXI(N*t)) 

S b A L * Y I f N ) ♦ ( 1»0 *aL)*YI f N • l ) 

P2 b 5+ xK*(Pl*$) 

350 Fl a ABS(PUS) 

E2 a ArS(P2»S) 

IFfEl*E2)400.400,4iO 
400 Y b 8 

GO TO 900 



«10 BT « (E1*AL)/(E1*ALt( 1 ,0"AL)*E2) 

Y a BT*P2+ C 1 « 0*BT ) *Pl 
POO ANS(IE)bY 
PIO continue 
RETURN 
END 
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SUBROUTINE 8TR A I T (F IN, X IN, NIN, XOUT , POUT, NQUT, NR I TE) 

C 

c this is a subroutine to determine the function fout at the 

C XCHJT iOCATTONS USING LINEAR INTERPOLATION from THE 
e < XIN, F IN5 table, 

c 

C *IN MUST EITHER BE IN aSCENDInG OR DESCENDING ORDER, 

^ C XOUT m a y be in any NUMERICAL ORDER, 

01 c NiN MAY BE 1 , IN N H I e m CASE TmE VALUE FjN ( ! ) IS ASSIGNED TO fFrjUT), 
C 

c***** this version is partially oourle precision ***** 

c 

DOUBLE precision xout,x 
real UNIT 

DIMENSION XIN(NIN) , FIN (NIN) ,XOUTCNOUT) ,FOUT(NOUT) 

. L I NIT (X, XI, X2,Fl,F2)«FlT(F2»Fl) ? ( X2-X 1) * ( X-X l ) 

c 

c 

c 

IF (NIN ,EQ, 1 ) GO TQ 600 
SIGNbI', 

IF(XIN(1) ,GT t XIN(2))SlGN3.1 , 

NUPRbN I N» 1 
00 500 NOb 1 , WOUT 

I^(SIGN*XOUT(NO) i GT l SIr,N*XIN(2) ) r,0 TO 100 


o 

ro 



FOUT(NO) s UNIT(XOUTtNO), X IN ( t ) , X IN ( 2 ) r F IN ( 1 ) , F I N ( 2 )) 

CO Tp 500 

100 IFf8IGN*X0UT(NO) t LT.SXpN*XXM(NUPR)) 50 TO 200 

FOOT f NO) * LINI T(XOUT( NO) #XIN(NUPR) »XIN( WIN) r FIN (NUP«) »FIM(NIN) ) 
50 TO 500 
200 CONTINUE 

00 300 NI«5,NUPR 

lF{SIGN*XOUTfNO) f G'F,8|GN*XlN(Nl)) GO TO 300 
FOUT(NO) b UINIT?X0UT(N0)^XIN(Ni-1) # XIN(N|), 
lFlN(NI-l)#FlNtND) 

GO TO 500 
300 CONTINUE 
C 

*ritf(nritf, n x i n 

X F0RM 4 T(jHi. IT 4 BLE NOT IN ASCENDING OR DESCENDING ORDER INI 
t r 1 SUB, STRAlTt/UXf IPf20,5)) 

STOP 13 
C 

500 CONTINUE 
RETURN 

600 DO 6lo NO® 1 ,W0UT 
610 rOUT(NQ)*FlN(!) 

return 

• **;• end ; 
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subroutine jntgrt 

(sJJl 

JJP. 

JJMaX 

# eta 

♦ stheta, 

F , ND 

* x, 

it w» 




lOUTR, NUFPT, V*UU) 

ccccccccc.ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccecc 
c 

c PROGRAM to INTEGRATE F (X) *SIN f X) FROM 0 TO PI by a trapezoidal RULE 
C MODIFIED To ALLOW FOR SPECIFIED DISCONTINUITIES, 
c THE DISCONTINUITIES MaY BE DISCONTINUITIES IN SLOPE OR VaLUE*. 


o o n n 



o n n n n n n n n r» r* ri c* o r* c> ri ni r* r* ci r* ? 


C IF ANY DISCONTINUITY Is a DISCONTINUITY IN VALUE AND THE LOCATION 

r OF THE OJScONIUITY COINCIDES WITH AN INTEGRATION POINT, THE PROGRAM 

C ASSUMES that the CORRESPONDING ELEMENT OF F IS FILLED RITH the LIMIT 

c average value, 
c 

C THE PROGRAM IS WRITTEN TO BE MACHINE INDEPENDENT 

C 

r 

JJ s NUMBER OF STATIONS TO USE IN THE INTEGRATION 

J JR e NUMBER OF VALUES OF F 

JJMaX * NUMBER OF VALUES OF ETa AND SThETa 

ETA s COS(ThETA) AT JJMAX STATIONS 

STHETA a S IN ( THETA) AT JJMaX STATIONS 

F e THE FUNCTION < VALUES AT JJR STATIONS 

ND s THE NUMBER OF DISCONT InUI T IES 

X s A VECTOR CONTAINING THE DISCONTINUITIES 

Y s THE X VECTOR AFTER REORDERING FROM HIGH TO LOW and . 

ELIMINATING VALUES OUT OF Range, DUPLICATED VALUES, 
and values WITHIN EPS OF ANY INTEGRATION STATION, 

OUTR e AN INTEGER GOVERNING THE OUTPUT OF THIS PROGRAM, 

NIEFT x THE NUMBER OF Y<*V ALUES , 

VALU * THE OUTPUT VALUE OF THE INTEGRAL 


■. all the input variables are returned 

the X-VALUES are completely unrestricted, but THE CORRECTIONS will 
NOT NECESSARILY BE DONE CORRECTLY in all cases, 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C,.,, .SPECIFICATION STATEMENTS 

double PRECISION ETa, SUM, PI, TD , STHETa 
INTEGER OUTR, OUT, W 
LOGICAL CONT, CONTI 

DIMENSION ETA(JJMAX), SThETa? JJMAX), F(JJR), XfND), Y f NQ) 

DIMENSION FT(a), FUNfa) 

EQUIVALENCE (YTemP, YDUM) 


oooocir»or*r*r>ooooorir*oocioT'»or*r»rir*r*oo 



c 

DATA PI /!', U1592653589793 DO/ 

DATA EPS 
C 

JJ, JJR , AND JJMaX must LIME up PROPERLY, ip THEY pQ NOT THE 

C i « 1 1 « program hill WRITE A MESSAGE TO that EFFECT and halt EXECUTION*. 

c 

IF <MOO(JJMaX*1, J JR* 1 5 ,Nfc, 0) GO TO 910 
IF<mod(JJR M* JJ M) .HE, 0) GO TO 9i0 
JJ1»JJ*1 

JR 1 «(JJR*1)/JJJ 
JR2 •CJJMAX + n/JJJ 
DELTA b PI/DFLOATCJJl) 

SUM»c 
OUT* OUTR 
C 

C • t « • , THE 50 LOOP DOES THE UNcORRECTED INTEGRATION 
C 

DO 50 J»1,JJ 
J2*J*JR2 

SUMsSUM+STHETAf J2)*FCJ1) 

50 CONTINUE 

VALU ■ DElTAbSUM 

THERE DISCONTINUITIES, TRANSFER TO 200. 

C 

IF(NO.OT,0) GO TO 200 
60 NLEPT»0 

IF (OUT t LT , 1 ) RETURN 
WRITE fW,i)VALU » JJ 
IF (OUT ,LT, ?.) RETURN 
IF(ND t £Q. 05 GO TO ?0 
WRITE (W, 23 
WRI.TE(W,5) X 

70 IF (OUT ,LT | 5) RETURN 

WRITE(W»<j) 

00 75 J«1,JJ 
J1 ■ J* JR l 



75 


J2aJ*JR2 

WRITE (W,5) ETA ( J2) , F(J1) 

RETURN 

c . .......... —w... ........ 

200 NtEpT * NO 

C 

c ,,,, .prior to correcting for discontinuities in slope and value the 

e EXTRANEOUS POINTS aRE ELIMINATED from CONSIDERATION, A TRANSFER 

c back to 60 will occur if there are no real discontinuities, 

c 

c 

C ARRANGING x in y in DESCENDING order 

c 

DO 210 Nb1,ND 
210 V ( N ) ■ X ( N ) 

DO 220 Nb 1 , NO 
DO 220 JaN, ND 
IE(Y(Ni p GE,Y(J))GO TO 220 
YTEMPbY C J) ‘ 

YtJ)«Y(N) 

Y(N) bYTEMP 
220 CONTINUE 

IF (NLEFT , EQ, 1 ) GO TO 300 
LMjNbI 
LMAx*ND-1 
C 

c,,.,, repeated values are eliminated next 
c 

230 DO 260 L«LMIN,LMAX 

IF(ABS(Y(L)-Y(L+i)).LT.EPS)GO to 250 
240 CONTINUE 
GO TO 260 

250 CALL CRUNCH (Y»L#LHAX) 

IF(L.GT,LMAX) fi 0 TO 260 

LMJNbL 

GO TO 230 

260 NLEFT 8 LMAX +l 
300 CONTINUE 
C 



C » * * * • TOO SMALL Y*-VAlUES ARE ELIMINATED NEXT 
C 

YDUMb EPS* 1 » 

DO 5 1 0 Iml'Hltn 
IFC Y(l 3 .LT.YOUMi GO TO 320 
310 CONTINUE 
GO TO 330 
320 NLfPTBL-i 

IF(NLEFT' # E9.0) GO TO 60 
330 IF(Y(NLEFT).GT,»YDUM5 GO TO 60 
400 IF (NLEFT. FQ. 1) Go To 500 
LMAXbNlEPT.1 

YDUMb-YDUM 

C l#lii TOO LARGE Y.VALUES are ELIMINATED next 

c 

410 IF CYCD.LE.YOUM) GO TO 440 
Call CRUNCHCY, 1,LMaX> 

GO to «10 

440 NLFFT B LMAXaI 
500 LMINbi 

y. Values nearly coincident hith integration stations are eliminated next 

e*"’ 

00 530 JbJR2, JJMAX.4R2 
£S B ETA(JJ 
DO 510 LbLMIN, NLEFT 
IF{ABSCY(U-I9ULT, EPS) GO To 520 
510 CONTINUE 
GO To 53o 

520 call CRUNCH (Y,L, NLEFT) 

IF(NLEfT,EG,0)GO TO 60 
1F(L,GT,NLEFT)G0 TO 540 

LMIN B L 
530 CONTINUE 

c . ..aT THIS POINT all the EXTRANEOUS DISCONTINUITIES HAVE BEEN 
ci:,,. ELIMINATED AND THERE ARE STILL SOME REMAINING, THE REMAINING DIS- 
C,. ,,, CONTINUITIES ARE STORED in (Y), 



n n n o non 


c 

540 KQUNT s 0 

CONT * , FALSE, 
V*LU!«VAUU 

IFCJJ.LT, 4) CO TO 800 


.♦t.fAT THIS POINT THE PROGRAM WILL MAKE THE NECESSARY CORRECTIONS 
•••••51 NOT BE ABLE TO DO THIS RIGHT UNDER CFHTaIN CONDITIONS. 

•■«»•!? t h ere are discontinuities in t m e first or last two intervals, 
MORE THAN one discontinuity in an interval, or discontinuities 
•••••I* ADJACENT INTERVALS, the corrections WILL not BE done right. 


FT(n»ETA(JR2) 

ET(2)«ETA(?*JR2) 

ETB)*eTA(5*jr2) 

CORRsO 

CONTI., false. 

JMAXs J J*3 
LMTN.J 

t, . 06 700 J.l'JMAX 

INDEV«(J*3)*JR2 
; ETt4)»ETA(INi5EX) 

DO 600 L* LMIN,NLEFT 
IF ( Y (L) ,GT,fT (2) ) GO TO 60Q 
IE(Y(L) .GT.eTB)) GO TO 610 
CONTI*'. FALSE. 

CO TO 6*0 
600 CONTINUE 
GO TO 79o 
610 INOEXaft 
J R 3 * J ♦ 3 

00 620 K*J,JR3 
K 1 #K* JR 1 
K2*K*JR2 
Ind.EX = INoEX + 1 

620 EUN(IN0Ex)*EfKl)*STHETA(K2) 
TO s DARCOS(DBLE(Y(L) ) ) 

02* TD/DELTA-oFLOAT( J+J) 

D3 si, - 02 



CORK ■ CnRR ♦ D2**2*(FUN(2)-FUN( 1 ) ) 
1 ♦ (2,*D2 • 1, )*(F!!Nt2)*FUM(3) ) 

COM * COM ,OR ,CONT 1 


CONTI* ; TRUE . 

KOUNTskQUNT* 1 

lF(kOUNT,EQ.NLEFT)GO TO ?90 

690 continue 

ET ( 1 ) a ET (2) 

FT(2)*ET{3) 

FT (35sfT(4) 

LMlNal 

700 CONTINUE 

790 VAIU*VALU +C0RR*DEtT A/2, 

*..—*#-* 

C,.,..ThIS is THE END OF THE COMPUTATION’. 
C»*t«»IS OUTPUT, 


c 

800 IF (OUT ,LT, 05 RETURN 

IF (KOUNT,NE,NUEFT)G0 TO » 1 0 
IF(CONT). GO TO 010 
s IFCJJ.lt. 9) SO TO 810 

IF (OUT.CT.O) SO TO 820 
return 

810 OUT ■MAXo (OUT# 2) 

820 HRlTf (w.fe) VALUl.VALU* JJ 
IF (OUT.IT ,2) RETURN 
IF (KOUNT.EQ.NlEFTl SO TO 830 
KOUNT * NLEFT • KOUNT 
WRITEOOT) KOUNT 
830 IF (.NOT'.CONT) SO TO 840 
HRITf(w,8) 

840 IF(JJ,GT,3)S0 TO 850 
WRITE (W, 9) 

850 WKITE(W,io) 

WRITg (W» 3 J X 

wrItf ( w # i ij 

WRITE (W » 3) (V(N),N»1,NLEFT5 
GO TO 70 


+ D3**2* (FUn(3)-FUN(4) ) 


THE REMAINDER OF THE PROGRAM 



o 

UJ 


910 WRITE (W, 12) JJ, JJR, JJMAX 
STOP 

I F0Rm a T(25H0VaLUE OF THE INTEGRAL IS ,F15,7 , 10X, 

HOh THE NUMBER Of INTEGRATION points USED * » |S) 

? PCRMAT(54H(5N0NE OF the x.values were CONSIDERED DISCONTINUITIES /• 
1 1 Th the X VECTOR is 1 
3 FORMAT ( 6F20,fe 1 

a FORMAT ( 1 OMO ETA , lflX, 1HF) 

5 FORMAT ( F12.6, FlJ,6) 

fc FORMAT ( 

IRIhOThE UNCORHECTED VALUE Of THE INTEGRAL a , F 1 5 » 7/ 

241 H THE CORRECTED VALUE Op THE INTEGRAL * >f!5,7/ 

341 H THE NUMBER Op INTEGRATION POINTS WAS a ,j7 ) 

7 F0RMaT(25HC- — WARNING — THERE WERE , 15, 

152h DISCONTINUITIES Which COULD NOT be CORRECTED FOR, ) 

8 FORMAT ( SOHO — WARNING — THERE WERE DISCONTINUITIES IN CONTIGUOUS- 

l INTEGRATION INTERVALS ) 

9 FORMAT (100H0— WARNING — NO DISCONTINUITIES aRE CORRECTED FOR WHE- 

IN There ARE LESS Than a INTEGRATION POINTS ) 

10 FORMAT (26H0O»IGINaL DISCONTINUITIES j 

II FORMAT (25H 0 REDUCED DISCONTINUITIES ) 

12 FORMAT (53H0JJ, JJR, A Nr> JJMAX ARE NOT PHOPpRLY RELATED, /- 

19h0 JJ b , 17/ 

19h JJR s , 17/ 

1 9M JJMAX ® 4 17/ 

13«HOEXECUTlON TERMINATED in INTgRT ) 

end 


END 

END 

end 

end 

end 

END 

END 

end 

end 

end 

end 

ENO 

END 

end 

end 

END 

END 

end 

end 

end 

end 

end 

end 

end 

end 

end 

END 

END 

end 

end 

end 

end 

end 

End 

end 

end 

end 

end 

end 

End 

end 

end 

END 

end 

END 

end 

end 

end 

END 

end 

end 

end 

END 

End 

end 

end 

end 

end 

end 

END 

end 

end 

END 

end 


'..RJStaRT RHt05 

SUBROUTINE CRUNCH (Y, L, LMaX) 
DIMENSION V C 1 ) 

DO 10 k * It LP*X 



10 


Y(K) s Y C K ♦ 1 ) 
LMAX ■ LMAX . 1 

return 

F Nf 5 


END 

END 

end 

end 

END 

End 

END 

END 

end 

END 

END 

end 

end 

end 

end 

END 

END 

END 

end 

end 

end 

END 

END 

end 

end 

END 

END 

end 

end 

END 

end 

END 

END 

END 

END 

end 

end 

END 

end 

end 

end 

END 

end 

end 

end 

END 

end 

END 

END 

END 

end 

END 

END 

end 

END 

end 

end 

END 

end 

end 

end 

end 

end 

end 


FUNCTION C 2 I(ETA,EtA 2 ) 

C # « » , THE FIRST 3 DERIVATIVES w.R.T. eta Vanish at eta2', 

AE 2 .ABSCETA 2 ) 

AEsABS(ETA) 

C2Ib1, 

IFeAE,LE.AE 2 jRETURN 

ZBfAE-AEZJ/Cl.-AEi) 

Z*Z*Z 

C2I«SQRT(l,-z)*n .+Z/2.) 

RETURN 

END 


END 

end 

end 

end 

end 

end 

end 

end 

end 

end 

end 

END 

END 

end 

END 

END 

END 

end 

end 

End 

end 

END 

end 

end 

end 

end 

END 

END 

END 

end 

END 

END 

end 

end 

end 

end : 

end 

END 

end 

end 

end 

end 

END 

End 

end 

END 

END 

end 

END 

end 

end 

end 

end 

end 

end 

end 

end 

end 

end 

end 

end 

end 

END 

end 


function Q ( X » Y ) 

C 

THIS FUNCTION IS FOR COMPUTING ACCURATELY THE V 4 LUE OF THE EXPRES* 

C« • • • • SI ON SGRT(X**?-fY**2} ♦ Y EVEN NHEN Y ,LT. 0, 4 ND A0S(X/Y) IS 

Ct.tt.MUCH smaller than i, 

c 

SIGN s tl, 

GO TO 5 



u UU tft U O <J 


ENTRY QhXNU5(X»Y) 

ENTRY HERE COMPUTES SQRTfY**2 » X**2) + V FOR THE SAME CONDITIONS, 
SIGN • 

IF { Y , GE , 0 , ) GO TO 10 
IF(ARS(X/Y) ,GT i ,15) GO TO 10 
Rb (X/y ) **2*$lGN 

TaYi. 0R»S EXPANSION OF -Y*(SQRTCl+R**2)-n, 

Qc -Y*R/2,*(l, " R/R,*(l, • R/2,*(l, * ,625*R))) 

RETURN 

10 Q a SORT f Y*Y ♦ SlGN*X*X) A Y 
RETURN 

End 


END 

end 

END 

end 

end 

end 

end 

END 

end 

end 

end 

end 

END 

END 

end 

end 

end 

end 

end 

end 

end 

END 

end 

END 

end 

end 

end 

end 

END 

end 

end 

end 

END 

END 

end 

END 

end 

End 

end 

ENp 

END 

END 

End 

end 

END 

end 

end 

END 

END 

END 

end 

END 

end 

End 

end 

end 

END 

END 

end 

END 

END 

END 

end 

end 


j- *8888888888888888888888888888888888888888888886888888888888888888688888888888 
o 


«ji ***** ***** 

***** DELETE CALLS TO cVRT,QREY,BcFIL ( and GEMFIL FOR OTMfR than ***** 

***** ames* version of the program ***** 

***** ***** 





